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ABSTRACT 

Baladi,  George  Youssef,  Ph.D.,  Purdue  University,  January  1968. 
Distribution  of  Stresses  and  Displacements  Within  and  Under  Long  Elastic 
and  Viscoelastic  Embankments.  Major  Professor:  Dr.  William  H.  Perloff . 

The  research  was  initiated  to  predict  stresses  and  displacements 
within  and  under  long  elastic  and  viscoelastic  embankments.  In  both 
cases,  the  embankment  was  assumed  to  be  homogeneous,  isotropic  and  con- 
tinuous with  the  underlying  material. 

As  the  first  step  in  the  development  of  a  suitable  theory,  the 
total  stress  was  defined  to  be  the  results  of: 

1  -  The  stress  due  to  the  geometric  configuration  of  the  boundary 

of  the  embankment. 

2  -  The  stress  due  to  the  weight  of  the  material. 

This  definition  was  used  together  with  the  Airy  stress  function  to 
represent  the  stresses  in  terms  of  analytic  functions. 

The  Schwarz-Christoffel  transformation  was  used  to  transform  both 
the  geometric  configureation  and  the  boundary  conditions  for  the  purpose 
of  matching  the  boundary  of  the  embankment,  with  the  stress  functions. 
The  Cauchy  integral  formula  was  applied  to  the  transformed  boundary  con- 
ditions to  yield  the  stresses  and  the  displacements. 

The  material  of  the  embankment  and  the  foundation  was  assumed  to 

behave  as  a- linear  viscoelastic  body,  the  time  dependent  solution  for 

the  stresses  and  the  displacements  was  obtained  by  applying  the  quasi- 

elastic  approximation,  proposed  by  Schapery  [32],  to  the  associated 

*  The  numbers  in  brackets  refer  to  references  given  at  the  end  of  the 
text. 
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elastic  solution  with  creep  functions  in  place  of  elastic  constants. 

Tine-dependent  material  properties  were  selected  from  the  results 
of  triaxial  creep  tests  on  a  compacted  clay  soil. 

The  numerical  results  of  both  the  elastic  and  the  viscoelastic 
solution  were  obtained  by  the  aid  of  an  U3I1  709^*  digital  computer  using 
a  Fortran  four  source  program. 

The  method  developed  permits  solution  of  a  wide  variety  of  plane 
strain  problems  with  variety  of  boundary  conditions. 

The  distribution  of  vertical  stress  underneath  the  elastic  embank- 
ment is  more  nearly  uniform  than  is  usually  assumed. 

The  horizontal  normal  stresses  was  found  to  be  significantly  lower 
than  those  computed  by  the  normal  loading  approximation. 

The  horizontal  shear  stress  acted  along  the  base  of  the  embankment, 
was  found  to  have  a  significant  value,  while  it  is  assumed  zero  by  the 
Eoussinesq  analysis. 

The  net  increase  in  vertical  normal  stress  within  the  compensated 
elastic  foundation  was  found  to  be  smaller  than  that  usually  assumed. 

The  viscoelastic  analysis  of  the  embankment  indicates  that,  for  the 
soil  investigated,  sixty  per  cent  of  the  displacements  which  occur  within 
the  first  ten  years,  is  due  to  the  elastic  effect  and  ninety  per  cent  of 
these  displacements  occur  in  the  first  month. 


INTRODUCTION 

The  distribution  of  stresses  and  displacements  within  and  under 
earth  embankments,  due  to  the  embankment  weight,  is  of  great  interest 
to  civil  engineers  in  a  variety  of  applications.  Consideration  of  de- 
formations within  embankments,  consolidation  of  underlying  compressible 
materials  and  other  engineering  applications  require  determination  of 
the  distribution  of  these  stresses  and/or  displacements.  Such  a  prob- 
lem has  been  under  active  investigation  by  various  investigators.  These 
investigations  are  mostly  concentrated  on  the  use  of  the  linear  theory 
of  elasticity.  However,  almost  all  these  studies  are  primarily  attempts 
of  approaching  a  solution. 

The  stress  distribution  within  a  body  depends  on  many  parameters, 
such  as : 

1  -  The  physical  properties  of  the  body,  which  may  be  a  function 

of  load  and  deformation,  and  the  stress  history. 

2  -  The  way  in  which  the  body  is  loaded. 

3  -  The  shape  of  the  loaded  area. 

k   -  The  geometric  configuration  of  the  body. 

Any  attempt  of  an  exact  solution  to  the  problem  must  consider  the 
above  parameters  properly. 

The  specific  objective  of  this  investigation  is  to  develop  a 
closed  form  solution  to  the  embankment  problem  to  predict  the  stress 
distributions  as  well  as  the  elastic  and  viscoelastic  displacements. 


It  is  assumed  herein  that  the  embankment  and  the  foundation  material 
with  which  it  is  continuous,  are  composed  of  homogeneous,  isotropic, 
linear  viscoelastic  material.  Further,  the  embankment  is  assumed  to 
be  sufficiently  long  so  that  plane  strain  conditions  apply. 

The  solution  is  obtained  by  transforming  the  region  of  the  embank- 
ment where  the  solution  is  unknown,  into  a  half-space  where  the  solution 
can  be  found.  Application  of  the  Cauchy  integral  formulas  to  the  boun- 
dary conditions  permits  determination  of  the  stresses. 


REVIEW  OF  LITERATURE 

Application  of  the  Theory  of  Elasticity  to  Soil  Media 
The  determination  of  stresses  and  displacements  within  a  loaded 
body  requires  assumptions  about  the  nature  of  the  loaded  material  and 
the  magnitude  of  the  resulting  displacement.  One  theory  available  for 
the  determination  of  such  stresses  and  displacements  assumes  that  the 
material  is  continuous,  homogeneous  and  isotropic  and  that  it  recovers 
its  original  shape  completely  when  the  applied  loads  are  removed.  The 
theory  based  on  these  assumptions  is  called  the  theory  of  elasticity. 
During  the  past,  the  theory  of  elasticity  has  found  considerable 
application  in  the  solution  of  engineering  problems.  But  application 
of  this  theory  in  its  general  form,  is  frequently  difficult,  because 
the  relation  between  stress  and  strain  is  generally  a  mathematical  equa- 
tion which  may  be  described  in  the  form  of  a  power  series.   Since  the 
behavior  described  by  such  a  series  is  nonlinear,  further  formulation 
presents  some  difficulties,  therefore  further  assumptions  are  intro- 
duced by  neglecting  all  the  terms  of  the  series  of  order  higher  than 
unity.   In  this  case,  the  strain  is  directly  proportional  to  the  ap- 
plied stress  and  can  be  described  by  Hooke's  law.  The  resulting  theory 
is  called  the  theory  of  linear  elasticity.  The  manner  in  which  the 
boundary  forces  are  distributed  within  the  soil  mass  depends  on  two  sets 
of  conditions:  first,  the  way  in  which  the  load  is  placed  on  the  soil, 
including  the  shape  of  the  load  over  that  area;  second,  the  physical 


properties  of  the  elastic  medium. 

With  a  few  exceptions  there  is  no  construction  material  whose  real 
mechanical  properties  are  simple  enough  to  be  acceptable  as  a  basis  for 
theoretical  analysis.  Therefore,  simplifying  assumptions  have  been  made 
regarding  these  properties  which  make  it  possible  to  compute  the 
stresses.  These  assumptions  are  always  to  a  certain  extent  at  variance 
with  reality. 

Taylor  [33]  said  for  the  theory  of  elasticity  to  be  applicable,  the 
real  requirement  is  not  that  the  material  be  necessarily  elastic,  but 
that  there  must  be  constant  ratios  between  the  stresses  and  strains. 
Therefore,  he  suggested  that  the  elastic  theory  may  be  applied  to  a 
non-elastic  soil  mass,  to  express  the  solution  in  mathematical  form, 
provided  such  constant  ratio  exist. 

Leonards  and  Narain  [2k]   conducted  bending  and  compressive  tests 
on  a  variety  of  soils,  both  undisturbed  and  remolded.  Figure  1  shows 
the  results  of  bending  stress-strain  curves  tests  performed  on  beams 
(3  in.  wide,  2  3/k   in.  deep,  and  22  1/8  in.  long)  of  compacted  limestone 
clay  from  Southern  Indiana.  Figure  2  and  Figure  3  show  the  results  of 
unconfined  compressive  stress-strain  curves  tests.  These  tests  were 
performed  on  cylindrical  specimens,  (l.k   in.  in  diameter  by  2.8  in.  in 
height)  on  the  soils  mentioned  before.  Figure  k   shows  stress-strain 
curves  for  compression  tests  on  both  undisturbed  and  recompacted  clay 
samples.  Dash  [10]  has  conducted  direct  tension  tests  on  consolidated 
clay  samples  (2.8  in.  long,  l.k   in.  in  diameter).  Figure  5  shows  the 
results  of  this  test.  These  results  of  bending  compressive  and  tensile 
tests  indicate  that  for  small  strains  the  resulting  stress-strain  curves 
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are  very  nearly  linear.  On  the  other  hand,  when  the  shear  failure  is 
imminent,  the  strains  are  far  from  being  proportional  to  the  stresses. 

An  indication  of  whether  theory  may  be  applied  in  a  given  situa- 
tion may  be  obtained  by  subjecting  soil  specimens  in  the  laboratory  to 
anticipated  field  stresses  and  plotting  the  elastic  effect  of  stresses 
against  strains.  If  straight  line  plots  are  obtained  the  theory  of 
elasticity  may  be  applied.  The  effect  of  time  may  also  be  important, 
viz.  Figures  2  and  3.  This  is  discussed  below. 

The  application  of  the  theory  of  elasticity  to  soils  where  the 
stress-strain  ratios  are  functions  of  depth  is  awkward  and  difficult. 
Holl  [18]  solved  the  problem  of  a  concentrated  force  acting  on  a  sur- 
face of  semi-infinite  mass  where  the  stress-strain  ratio  varies  in  the 
form: 


E  =  E  zX 


where  E  is  the  modulus  of  elasticity  at  z  =  1  and 


A  +  z  =  i 

where  u  is  Poisson's  ratio. 

However,  Taylor  [33]  suggests  that  in  many  soils  this  variation  is  not 

so  large,  and  that  the  overall  action  is  about  the  same  as  in  the  hypo- 
thetical case  in  which  the  stress-strain  relationships  are  constant  and 
equal  to  the  average  value. 

The  theories  of  behavior  of  soils  under  stress,  are  very  compli- 
cated and  the  geometric  configuration  of  the  boundary  of  the  soil  me- 
chanics problems  are  not  simple,  therefore,  the  use  of  the  linear  theory 
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of  elasticity  has  been  limited  to  problems  which  may  be  idealized  to 
semi-infinite  solids  with  loads  on  their  boundary  planes.  Terzaghi  [3*+] 
has  pointed  out  that  in  applying  the  theory  of  elasticity,  the  investi- 
gator needs  to  know  (a)  how  to  compute  the  physical  properties  of  the 
soils;  (b)  how  the  results  thus  obtained  compare  with  experience.  Never- 
theless, the  theoretical  solution  can  give  us  insight  into  the  problem, 
clarify  both  the  limit  of  the  solution  and  the  boundary  conditions;  and 
the  stress  components  at  points  in  the  region  can  be  computed  from  known 
equations  defining  the  stresses. 

Viscoelastic  Material  Representation 
A  viscous  material  is  one  whose  constitutive  relations  are  time- 
dependent.  A  viscoelastic  material  is  one  which  combines  elastic  and 
viscous  response.  As  in  elasticity,  the  materials  are  usually  assumed 
to  be  linearly  viscoelastic  in  order  to  simplify  their  representation 
and  to  be  able  to  apply  superposition. 

Whether  or  not  a  real  material  can  be  described  as  linearly  visco- 
elastic must  be  determined  experimental ] y  [21,  38]. 

For  certain  problems  in  linear  viscoelasticity  the  correspondence 
principle  can  be  used  to  calculate  a  time-dependent  viscoelastic  re- 
sponse from  the  solution  of  an  "associated"  elastic  problem.  This  prin- 
ciple was  deduced  by  Lee  [22]  for  isotropic  media  and  later  by  Biot  [2] 
for  anisotropic  material.  The  basis  of  this  principle  is  that,  with 
zero  initial  conditions,  the  Laplace  or  Fourier  time-transformed  visco- 
elastic field  equations  and  boundary  conditions  are  formally  identical 
with  the  equations  for  an  elastic  body  of  the  same  geometry.  Thus 
transformed  solutions  can  be  obtained  by  standard  elastic  analysis, 
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and  then  inverted  to  obtain  the  time -dependent  response.  Lee  and  Roger 
[23]  and  Valanis  and  Lianis  [37]  have  suggested  some  numerical  methods 
which  do  not  involve  transforms,  and  are  applicable  when  the  relation 
between  the  stress  and  strain  is  expressed  in  the  form  of  convolution- 
type  integrals,  involving  experimentally  measurable  creep  or  relaxation 
functions.  The  combination  of  this  stress-strain  equation  with  the  re- 
maining equations  of  viscoelasticity,  leads,  in  general,  to  integro- 
differential  equations  governing  both  space  and  time  dependence  of 
stresses  and  displacements.  The  solutions  of  these  equations  are  ap- 
proximated by  using  finite-difference  integration  with  respect  to  time 
[37]. 

Schapery  [32]  has  proposed  a  method  which  makes  use  of  the  hered- 
itary form  of  the  stress-strain  equations.  It  consists  first  of  cal- 
culating the  response  of  a  fictitious  elastic  body  which  is  geometri- 
cally similar  to  the  actual  viscoelastic  body,  has  the  same  mechanical 
loading,  and  has  elastic  properties  which  are  numerically  equal  to  the 
instantaneous  relaxation  moduli  or  creep  compliances  associated  with 
the  viscoelastic  problem;  these  mechanical  properties  are  expressed  as 
functions  of  real  time  in  isothermal  problems.  Two  advantages  of  this 
method  are  that  for  many  problems  good  accuracy  is  achieved  with  only 
the  quasi-elastic  approximations,  the  difference  between  the  exact  vis- 
coelastic solution  and  this  method  can  be  neglected;  and  it  is  poten- 
tially applicable  to  nonlinear  as  well  as  linear  problems  [32], 

Review  of  Previous  Contributions  to  the  Specific  Problem  Considered 
Elastic  theory  has  been  used  since  1920  to  investigate  stresses 
and  deflection  under  long  embankments.  At-  that  time  Carothers  [7] 
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determined  the  stresses  within  an  homogeneous,  isotropic  elastic  half- 
space  resulting  from  a  "long  embankment"  loading.  It  was  assumed  that 
the  embankment  can  be  represented  by  a  series  of  independent  line  loads 
applied  normal  to  the  boundary  with  a  magnitude  which  varies  in  such  a 
manner  that  it  graphically  resembles  an  embankment  cross-section,  as 
shown  in  Figure  6(b).  He  applied  superposition  on  equations  developed 
by  Boussinesq  [k]   and  Flamant  [15].  These  results  were  presented  in 
tabular  form  by  Jurgenson  [19].  In  19^+0,  Newmark  [26,  27]  used  the 
same  assumption  above  to  develop  his  solution  which  was  superimposed 
later,  by  Osterberg  [29],  to  obtain  an  influence  chart  for  the  deter- 
mination of  the  magnitude  of  vertical  stresses  induced  in  an  elastic 
half-space  by  a  long  embankment  loading  with  a  variety  of  cross-sections. 
These  approaches  to  the  problem  do  not  consider  the  elastic  characteris- 
tics of  the  embankment  or  the  shearing  stresses  at  the  interface  between 
the  embankment  and  the  foundation  due  to  elastic  body  restraint.  Thus, 
these  solutions  do  not  describe  the  real  physical  situation. 

Trollope  [36],  and  Davis  and  Taylor  [11]  computed  the  distribution 
of  stresses  underneath  a  granular  embankment  resting  on  a  foundation 
which  yielded  an  arbitrary  amount.  No  attempt  was  made  to  compute  the 
amount  of  foundation  movement  which  would  be  created  by  the  embankment. 
The  impression  given  by  these  methods  is  that  the  nature  of  the  problem 
is  so  complex  that  a  large  number  of  arbitrary  assumptions  must  be  in- 
troduced. In  particular,  the  assumption  of  limiting  equilibrium  is  a 
very  discouraging  feature  of  this  type  of  approach  to  the  problem. 

Dingwall  and  Scrivner  [12]  used  a  finite  difference  method,  to 
present  an  approach  to  an  elastic  embankment  resting  on  a  rock 
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(a)     Long     Symmetric     Elastic     Embankment 
Continuous  With     Foundation 
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(b)    Normal    Loading    Approximation 


Figure  6.       Problem    Considered 
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foundation.  This  method  was  used  by  Carlton  to  study  an  elastic  em- 
bankment continuous  with  an  elastic  foundation.  Zienkiewicz  [40,  43] 
applied  a  finite  difference  method  to  find  the  stress  distribution  with- 
in a  gravity  dam  resting  on  an  elastic  foundation  of  different  modulus. 

Zienkiewicz  and  Cheung  [Ul,  k2]   applied  the  finite  element  method 
to  find  the  stresses  within  buttress  dams  resting  on  an  elastic  founda- 
tion. This  method  has  been  used  by  Clough  and  Chopra  [9]  and  Finn  [13] 
to  the  study  of  a  triangular  dam  on  a  rigid  foundation  and  a  rock  slope 
continuous  with  its  elastic  foundation,  respectively.  Brown,  [6]  and 
Goodman  and  Brown  [17]  investigated  the  case  of  a  long  elastic  slope 
constructed  incrementally;  it  is  not  clear  to  what  degree  their  results 
are  influenced  by  the  fact  that  the  compatibility  equations  are  not 
satisfied  by  their  solution  method. 

All  the  solutions  obtained  by  applying  numerical  methods  are  ap- 
proximate and  were  either  restricted  to  a  single  embankment  cross- 
section  or  a  complete  stress  picture  was  not  obtained. 

In  connection  with  stability  analysis,  Brahtz  [5]  developed  a  meth- 
od for  analyzing  the  stability  of  earth  dams.  Brahtz  made  use  of  equi- 
librium equations  in  addition  to  the  assumption  that,  on  the  line  of 
symmetry  of  the  embankment,  the  vertical  and  horizontal  normal  stresses 
were  proportional.  Brahtz  assumed  that  the  boundary  was  free  from  ex- 
ternal and  internal  stress.  The  stress  distribution  results  from  this 
method  appear  to  be  arbitrary,  because:   (a),  the  only  basis  of  the 
suggested  formula  is  static  equilibrium,  which  is  a  necessary,  but  in- 
sufficient, requirement;   (b),  the  boundary  conditions  are  only  par- 
tially considered;  and  (c),  the  notion  of  constant  ratio  between  the 
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horizontal  and  vertical  normal  stresses  acting  at  every  point  of  the 
center  of  symmetry  is  not  properly  justified. 

Glover  [16]  has  considered  the  stability  of  an  earth  dam  with  tri- 
angular cross  section,  by  dividing  the  foundation  into  a  series  of  elas- 
tic and  plastic  regions.  The  alternating  plastic  and  elastic  zones  are 
assumed  bounded  by  straight  lines.  In  each  zone,  the  components  of 
stress  are  supposed  to  vary  linearly,  and  continuity  of  stress  is  as- 
sumed at  the  boundaries.  This  approach  introduces  a  large  number  of 
arbitrary  assumptions,  which  is  a  very  discouraging  feature  of  this  type 
of  approach  to  the  problem  of  stability  of  slopes. 

Terzaghi  [3^]  described  an  effort  by  Rendulic  [31]  to  evaluate  the 
shear  stresses  transmitted  to  the  foundation  material  by  an  embankment. 
The  shearing  stresses  acting  on  the  base  of  the  embankment  were  com- 
puted by  assuming  that  the  embankment  material  was  on  the  verge  of 
failure  and  thus  the  shearing  resistance  of  the  embankment  was  fully 
mobilized  in  order  to  maintain  equilibrium.  This  method  applies  the 
Mohr-Coulomb  theory  to  the  equations  of  equilibrium.  Consequently, 
Terzaghi  [3^]   suggested  that  the  magnitude  of  the  computed  shear  stress 
at  the  base  was  likely  to  be  lower  than  the  actual  in-situ  stresses. 

Finn  [lU]  suggested  the  use  of  the  Schwarz-Christoffel  transforma- 
tion to  map  the  embankment  region,  where  the  solution  is  very  difficult 
to  get,  into  a  semi-infinite  region,  where  the  solution  can  be  easily 
obtained.   However,  he  did  not  carry  out  the  suggested  procedure. 

Scope  of  the  Problem  Considered  Herein 
As  stated  above,  the  need  for  an  exact  solution  to  an  embankment 
with  its  foundation  is  evident.  In  this  investigation,  a  general  exact 
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method  of  solving  an  elastic  embankment  continuous  with  the  underlying 
material  is  developed.  In  fact,  it  is  believed  that  this  method  can  be 
applied  to  solve  many  plane-strain  problems,  under  statics  as  well  as 
dynamic  loads.  A  quasi-elastic  approximation  is  applied  to  obtain  the 
time-dependent  stresses  and  displacements  for  a  viscoelastic  material. 
The  material  properties  for  the  viscoelastic  response,  are  evaluated 
from  laboratory  test.  Numerical  results  and  influence  charts  for  a 
variety  of  embankment  cross-sections  are  presented.  Some  numerical  re- 
sults for  a  foundation  are  presented  to  show  the  potential  of  this 
method . 
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PART  I 
ELASTIC  ANALYSIS  OF  THE  EMBANKMENT 
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THE  EQUATIONS  FOR  THE  PLANE  STRAIN  THEORY  OF  ELASTICITY 

Introduction 

Let  u,  v  and  w  be  the  components  of  the  displacements  in  the  x,  y 
and  z  directions,  respectively.  A  body  is  said  to  be  in  a  state  of 
plane  strain,  parallel  to  the  plane  oxy,  Figure  7»  if  the  component  w 
is  equal  to  zero  and  if  the  components  u  and  v  are  independent  of  z. 
The  problem  of  the  determination  of  stresses  and  displacements  of  a  long 
elastic  embankment  continuous  with  its  foundation  can  be  formulated  an- 
alytically by  using  the  equations  of  plane  strain.  The  solution  will 
be  unique  if  it  satisfies  the  equilibrium,  the  compatibility  equations, 
and  the  boundary  conditions  [28]. 

The  equations  of  equilibrium  for  plane  strain  are: 

do     dt 
dx     or 

(1) 
bi  do 


in  which  pg  is  the  unit  weight  of  the  material,  o  and  o  are  the  total 

°  '  x      y 

normal  stress  components  parallel  to  the  x-  and  y-axes,  respectively, 
and  t   is  the  shear  stress  parallel  to  the  x-  and  y-axes. 

The  compatibility  equation  for  plane  strain  and  for  constant  unit 
weight  is : 
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Figure  7.        Two    Dimensional    Embankment 


Figure  8.         An    Element    from   the    Boundary 
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^(o  +  a  )  =  0  (2) 


where  v  is  the  harmonic  operator. 

The  boundary  stress  equations,  Figure  8,  are: 


X  =  a     cos(n,x)  +  x   cos(n,y) 

a  Ay 


(3) 


Y  =  t   cos(n,x)  +  or  cos(n,y) 


where  X  and  Y  are  the  x  and  y  components,  respectively,  of  a  distri- 
buted surface  force  per  unit  area,  cos(n,x)  =  ^  and  cos(n,y)  =  -  —  . 

Definition  of  the  Stresses 
Define  the  total  stress  within  the  region  under  consideration  as 
the  following: 

x     X     l-|i  " 


oy  =  o'   +  77  CO 


1         m   T » 

xy    xy 


where  a',  a'   and  x'  are  the  stresses  due  to  the  shape  of  the  region 

x'  y     xy  .  *  B 

only  and  u  is  Poisson's  ratio.  Therefore: 


o'  =  a'  «=  x'  =  0 

x    y    xy 


for  a  semi-infinite  region  with  no  boundary  loading.  Substituting 
equations  (h)   into  equations  (l)  and  (2)  gives: 
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S°x   ^Tx 
bx.         by 

(5) 

bo'       oV 

_JC  +  _^EZ=  o 

by        dx 


and 


which  are  the  equilibrium  and  compatibility  equations  for  a  weightless 
material . 

Stress  Function  for  Plane  Strain 
The  first  of  equations  (l)  represents  the  necessary  and  sufficient 
condition  for  the  existence  of  some  function  B(x,  y)  such  that  [35]: 

**&*■    -V   .   ^-&-^   (7) 

The  second  of  the  equations  (l)  is  the  necessary  and  sufficient 
condition  for  the  existence  of  some  function  A(x,  y)  such  that  [35]: 

M^ll    =  a       _  yy   ,      *(*,  7?  =  .T  (8) 

bx  y  by  xy 

Comparison  between  equations  (7)  and  (8)  for  t   shows  that: 

&(*>  y)  e  Mxt  y)  (9) 

by  bx. 

whence  follows  the  existence  of  some  function  U(x,  y)   such  that: 

A(x,   y)   .     M^ll  ,  B(x,   y)   =     M^JLl  (10) 
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Substituting  these  values  for  A(x,  y)  and  B(x,  y)  into  equations  (7) 
and  (8),  it  can  be  seen  that  the  stress  can  be  expressed  as: 


i&Si  &  ♦ 


_H_ 


^2       ^  ^ 


o 


y  ax2 


.  ,  ^(x,  y) 
xy       dxdy 


or,  in  terms  of  the  "geometric"  stress: 


^A  g^   +   ,y  (11) 


°x   "      .2 


°y  "    ,2 


^fejri  (12) 

dx2 


yfy*.  y) 


xy       axdy 

This  fact  was  first  noticed  by  G.  B.  Airy  (1862).  The  function 

U(x>  y)  is  called  a  stress  function,  or  the  Airy  stress  function  [35]. 

Assume  that  the  stress  components  a  ,   a    and  t   are  single-valued 

*  x'  y     xy       * 

and  continuous  functions  with  continuous  derivatives  up  to  and  includ- 
ing the  second  order  throughout  the  region  occupied  by  the  body.  There- 
fore, the  function  U(x,  y)  must  have  continuous  derivatives  up  to  and 
including  the  fourth  order  and  these  derivatives,  from  the  second  order 


2k 


onwards,  must  be  single-valued  functions  throughout  the  region  of  the 
embankment . 

If  U(x,  y)  has  these  properties,  the  function  a  ,  a  and  t  will 
satisfy  the  equilibrium  equations.  However,  for  the  uniqueness  of  the 
solution,  the  compatibility  equation  (2)  must  also  be  satisfied.  Rec- 
ognizing that: 


leads  to: 


or 


0x +  °y  -  ^X»  y) +  xk  ™ 


VUU(x,  y)  =  0  (13) 


*\*^  Y)     ,  c  ^Vx.  y)  t  aV«.  y)  ^  Q 

&x         Sx  Sy        oy 


This  equation  is  the  biharmonic  equation  and  its  solution  consists  of 
biharmonic  functions,  which  are  the  only  functions  that  satisfy  the  bi- 
harmonic equation,  the  derivatives  of  which  are  continuous  up  to  and 
including  the  fourth  order  and  are  single-valued  throughout  the  region 
under  consideration  [35]. 
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COMPLEX  REPRESENTATION  OF  THE  GENERAL  SOLUTION 

OF  THE  PLANE  STRAIN  PROBLEM 

INCLUDING  THE  WEIGHT  OF  THE  MATERIAL 


Introduction 
As  mentioned  above,  the  problem  under  consideration  is  formulated 
in  terms  of  a  biharmonic  differential  equation.  Direct  solution  of  such 
an  equation  is  generally  awkward  and  difficult.  Therefore,  the  trans- 
formation of  the  biharmonic  equation  to  a  harmonic  one  by  use  of  analyt- 
ic complex  function  is  expedient.  This  is  done  below  by  expressing 
stresses  and  displacements  as  complex  harmonic  functions  in  a  manner 
similar  to  that  developed  by  Muskhelishvili  [25]. 

Complex  Representation  of  the  Stress  Function 
Introducing  the  following  notation: 

sAj(x,  y)   -  O    +  a     -  1-^77  =  0'+  a'  =   P(x,  y)   (Ik) 


then 


V  U(x,  y)  =  v^P(x,  y)  =  0  (15) 


Thus  the  function  P(x,  y)  is  a  harmonic  function. 

Let  Q(x,  y)  be  a  harmonic  function,  conjugate  to  P(x,  y);  i.e.,  a 
function  satisfying  the  Cauchy-Riemann  conditions 


dP(*,  y)    .    aq(x,  y) 

cbc  dy 


(continued) 
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ap(x»  y)  =  _  ^(x;  y)  (l6) 

dy         dx 

This  function  is  determined  for  a  given  P(x,  y)  apart  from  an  arbi- 
trary constant  term.  Then  the  following  expression: 

W(z)  =  P(x,  y)  +  iQ(x,  y)  (17) 

will  represent  a  function  of  the  complex  variable  z  =  x  +  iy,  harmonic 
in  the  region  under  consideration. 
Furthermore,  put 

♦(z)  =  p(x,  y)  +  iq(x,  y)  -  J  /  W(z)dz  (l8) 

The  derivative  of  equation  (18)  with  respect  to  z  is: 

♦  •(Z)    =  MllJll  +   i  MZU±  =  J   [P(x,    y)  +iQ(x,    y)]    (19) 
dx  by 

Whence,  noting  that  by  the  Cauchy-Riemann  conditions: 

fo(x?  y)  =  ssl(x?  y)      »      dp(*>  y)  =  .  dq(*»  y)  (20) 

dx  cty  dy  cbc 


One  obtains : 


Wx>  y)  =  ^a(x;  y)  _  J  P(x  v\ 

dx  oy  5     k   '  Y> 

MemjA  =  .  M^ud  m  .  i  Q(x,  y) 


(21) 


dy  dx 

Thus: 
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Therefore : 


P(x,  y)  =  k  MLUd  s  k   **(*,  Y)   =  v2u(X)  y)    (22) 
ox         dy 


v2u(Xj  y)  .  2  fr(x,  y)  _  2  dq(x,  y)  =  Q  (g3) 

dx  dy 


This  leads  to: 

^[V(x,   y)  -  xp(x,  y)  -  yq(x,  y)]  =0  (2k) 

The  solution  of  this  equation  is  a  harmonic  function  in  the  region  under 
consideration.  Let  p, (x,  y)  be  such  a  harmonic  function,  then: 

U(x,  y)  =  xp(x,  y)  +  yq(x,  y)  +  P;L(x,  y)  (25) 

Now,  let  X(x,  y)  denote  the  function  of  the  complex  variable 
z  =  x  +  iy,  the  real  part  of  \vhich  is  p,  (x,  y)  and  the  imaginary  part 
of  which  is  a  function  q,  (x,  y),  which  is  the  harmonic  conjugate  to 
p1(x,  y) .   Substituting  this  into  equation  (25)  gives: 

U(x,  y)  =  xp(x,  y)  +  yq(x,  y)  +  ReX(z)  (26) 

Combining  xp(x,  y)  with  yq(x,  y),  gives: 

U(x,  y)  =  Re[(x-iy)p(x,y)  +  i(x-iy)q(x,y) ]  +  ReX(z) 
Substituting  x-iy  =  z  leads  to: 

U(x,  y)  =  Re|  z[p(x  ,y)  +  iq(x,  y)]  +  X(z)| 
Recognizing  that  p(x,  y)  +  iq(x,  y)  =  <t>(z);  gives: 
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U(x,  y)  =  Re[z  ♦(z)  +  X(z)]  (27) 

From  the  properties  of  the  complex  function  (Churchill  [8]),  one  can 
obtain: 


2U(x,  y)  «=  z  *(z)  +  z  *(z)  +  X(z)  +  X(z)  (28) 

The  partial  derivatives  of  this  equation  with  respect  to  x  and  y, 
respectively,  are : 

2*Uix>   y)  s  ♦(z)  +  7*«(z)  +  *00  +  z*»(z)  +  X'(z)  +  X'(z)   (29) 
Sx 

2^X>   y?  =  i[-*(z)  +  I*(z)  +  ♦(!)"  -  z*'(z)  +  X«(z)  -  X'(z)   (30) 

Adding  equations  (29)  and  (30),  after  multiplying  the  second  one  by  i, 
gives : 

M^il  +  i  M^il »  ,kz)  +  z^77y  +  xT^T  (31) 

dx  dy 

Let 

X'(z)   =     5g5l.t(s) 

Then 

Mis*_xl  +  i  MiS2_ii  =  ,(z)  +  zTmx)  +  tOO  (32) 

dx  by 

Differentiating  equation  (29)  with  respect  to  x  and  equation  (30)  with 
respect  to  y  and  adding,  leads  to: 

nAj(x,  y)  =  2[**(z)  +  ♦•(z)]  =  k   Re[*'(z)]  (33) 
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Equation  (33)  shows  that  vU(x,  y)  is  completely  determined  by  the 
real  part  of  the  function  *'(z).  It  is  also  convenient  at  this  point 
to  recognize  that: 

vAj(x,  y)  =  ax   +  ay  -  ^  yy   =  0^  +  cy  =  k   Re[*»(z)]  (3k) 

Complex  Representation  of  Stresses 
It  is  shown  in  Figure  8  that  the  force  (Xds,  Yds)  acting  on  an 
element  ds,  will  be  understood  to  be  the  force  exerted  on  the  side  of 
the  positive  normal.  Recalling  equations  (3) 


X  =  a     cos(n,x)  +  t   cos(n,y) 

A  Ay 


Y  =  Txy  cos(n»x)  +  ay  cos(n,y) 


(35) 


Combining  equations  (ll)  with  equations  (35),  leads  to: 


.2,,,    N  -.2, 


—   d  U(x,  y)    ,    v   c*  U(x,  y)    /    \    u        /    \ 
X  =  v  >  *  <    cos(n,x) v-  *  *  1    cos(n,y)  +  r~—  yy   cos(n,x) 

(36) 


ay2  ^x  *y 


Y  =  .  ^{x,   7)   cos(n,x)  +  ^U(x?  y?  cos(n,y)  4  yyC^y) 
^x  cV  ^2 


Where : 


:os(n,x)  >  cos(s,y)  m  ^ 


cos(n,y)  =  -cos(s,x)  =  -  — 


s  is  the  positive  direction  of  the  tangent, 
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Introducing  these  values  into  equations  (36),  gives 


1     ds 


or  in  complex  form: 


ds  L       oy  5x       J       ^  1>H  ds  dsj 

(X  +  iY)  =  -  id  pU(x>  ^  +  i  M*i_zl]  .  iy7  [ax  +  i  JL  dy] 
Substituting  equation  (32)   into  equation  (38)  gives: 


(38) 


(X  +  iY)ds  =  -  id 


U(z)   +  z*'(z)   +  ^(z^-iJdx  +  i  -H_  dyj    (39) 


First  let  the  element  ds  have  the  direction  of  the  axis  oy. 

Then 

ds  =  dy 

dz  =  idy 

dz  =  -  idy 

X  =  0 

X 

I   =  T 

xy 

Substituting  these  values  into  the  above  expression: 

°x  +  l   Txy  =   *'(z)    +   *'(z)    "  Z*"(z)    "  *'(z)    +  T^I     y        (U0) 
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Next,  let  ds  have  the  direction  of  the  axis  ox.  Then: 
ds  =  dx 
dz  =  dz  =  dx 

X  =  -  T 
Y  =  -  a 

y 

From  the  sane  expression  after  multiplying  both  sides  by  i: 


ay  -  i  t   =  *'(z)  +  ♦'(z)  +  zO"(z)  +  t|r'(z)  +  7y      (Ul) 


Equations  (UO)  and  (hi)   are  the  required  expressions  for  the  stress 
components.  Adding  and  subtracting  these  formulas  and  replacing  in  the 
latter  case  i  by  -i,  the  following  is  obtained: 


o 

X 


+  ay  =  2[*'(z)  +  ♦•(z)]  +  ^  yy 

=  k   Re|*(z)]  +  ~  7y 

=  2[0(z)  +  *(!)]  +  ■—  7y  (1+2) 


where : 

<J>(z)  =  ♦,(z) 
and: 

a     -  a     +   2ii  =   2[z4(z)  +  f(z)]  +  —H  /y  (U3) 

y    x      xy    l   v  /    i\  /j   2._(j  /Jr  v  -" 

where : 

f(z)  =  ^(z) 

Because  $(z)  and  $(z)  are  analytic  functions  in  the  region  under 
consideration  and  have  continuous  derivatives  up  to  and  including  the 
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second  order,  the  quantities 


a     +  a     -   = —  ry 

x    y   1-u  /J 


1-2 

O  -  0  +  2ix   -  •: yy 

y    x     xy   1-u 


are  also  analytic  functions  in  the  same  region  and  have  continuous  de- 
rivatives up  to  and  including  the  second  order. 

Complex  Representation  of  the  Displacements 
The  stress-strain  relationships  for  plane  strain  are: 


but 


~  e     =  a'   -  u(o'   +  a') 
Hu     x         xx  y 


E  ,  /    ,  ,\        l-2u 

-r— -  e     =  a'   -  p(a'   +  o1)   +  ■= — -  7y 
1+m    y        y  x        y'      l-u 


°x +  °v  =  ^u<x>  y)  =  p(x>  y) 

A        y 


(MO 


then 

.2 


E  ^U(x,   y)  ?..,  x 

1+H  £x  "  ^2  "   ^U(x>   y) 


(*»5) 


E  d2U(x,   y)  _2TI,  v    .    l-2|i 

^  ey =       L*       '  **    ™ 

where  E  is  the  elastic  modulus. 

Replacing  in  the  first  of  equations  (kk)       U^Xj!,  ?'   by  P(x,y)  -  "^(x?  y) 

5y  Sx 
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and  in  the  second  equation  *  '  *<   by  P(x,y) *  ■'  *■*■,   leads  to; 

bx  by 


2G  ex  =  -  ^x'/)  +  (m)  P(x,y) 
bx 

2G  e      ^%jd  +  (1-n)  P(x,y)  +  i^i  7y 
y       Sy2  ^ 

Where: 


G  =  p/vr  \   is  the  shear  modulus 


Substituting  the  value  of  P(x,  y)  from  equation  (21)  into  the  above 
equation: 

2Gcx  =  .  52u(x>  y?  +  m-^)  fr<x>  y) 

dx  dx 

(hi) 

2G  e   .  .  yW,  y?  +  ud.,)  **(*,  y)  +  JsSk 

y       Sy2  Sy      ^ 

Integrating  the  first  equation  with  respect  to  x  and  the  second  with 
respect  to  y: 


2Gu  =  -  5Ufo  r)   +  kd-v)  p(x,   y)   +  f(y) 


2 


2Gv  =  .  M2b_yl  +  m.^  q(x,  y)  +  ^h  zz_  +  g(x) 

by  ^      2 


To  determine  the  functions  f(y)  and  g(x),  differentiate  the  first 


3h 


equation  with  respect  to  y  and  the  second  one  with  respect  to  x.  Add- 
ing the  results,  gives: 

2Gj|i  +  |xj  .  .o  a2"(*,  y)  +  k(i-v),  ^(x?  y)  +  **(*>  y>  I 

L  y       c    -■  Sx  Sy  L       Sy  Sx       -• 

+  f (y)  +  g'(x)  (U9) 


Recognizing  that 


2Gi|Ugl|   =  .o     "^   ^  (50) 


and 


[ay  +  dxj   =  "  ^   Sy 


^(x?  y)  +  ^(x>  y)  =  0 

dy  dx 

results   in 

f (y)  +  g'(>0  «  o 

Because  f'(y)   is  a  function  of  y  only  and  g'(x)  is  a  function  of  x 
only,  this  cannot  be  possible  unless : 

f(y)  =  cy  +  c1 

(51) 

g(x)  =  -ex  +  c2 

where  c,  c,  and  c„  are  arbitrary  constants.  Substituting  these  values 
into  equations  (U8)  leads  to  expressions  for  u  and  v  which  are  definite 
apart  from  terms  of  the  form: 
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^  =  cy  +  cx       ,        \'1   =  -ex  +  c2  (52) 

These  terms  express  only  rigid  body  displacement  (in  the  oxy  plane) 
and  they  do  not  influence  stresses  and  strains  [25].  The  constants, 
c,  c.  and  c?  attain  definite  values,  if  one  assumes  as  given  the  values 
of  the  components  of  initial  rigid  body  displacements  and  rotation  of 
the  region  under  consideration. 

Therefore,  equations  (k8)   can  be  written: 


2Gu  =  -  M*2-ll  +  l4(l.^)p(x,   y) 
5x 

2Gv  .  .  M^ll  +  m_,)q(:,,   y)   +  "K  ZL- 
By  ^        2 


(53) 


also 


2G(u  +  iv)  = . ,  au(*»  y?  +  i  Maud,  +  u(i-M)*(2) 

L   Sx  Sy   -1 

Combining  equations  (32)  and  (5*0  gives: 


2G 


\(vl   +  iv)  =  (3-Un)4.(z)  -  z  ♦•(z)  -  \|r(z)  +  i  y5^  ^y2  (55) 
which  determine  the  displacement  in  terms  of  the  complex  functions 


4>(z)  and  \|r(z) , 
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APPLICATION  OF  CONFOKMAL  MAPPING 


Introduction 


In  order  to  find  the  value  of  the  stress  function  derived  above, 
the  boundary  of  the  problem  has  to  be  smooth.   Since  the  boundary  of  the 
embankment  is  composed  of  several  straight  lines  intersecting  at  angles 
between  zero  and  it  radian,  i.e.,  the  boundary  is  not  smooth,  then  it 
must  be  transformed  to  a  smooth  curve.  A  conformal  mapping  is  used  to 
obtain  a  closed  form  solution. 

The  Transformation  of  the  Geometric  Boundary  into  a  Half-Space 
The  boundary  of  the  embankment  is  a  polygon  with  one  of  its  ver- 
tices at  infinity.  Thus  the  Schwarz-Christoffel  transformation  will 
map  it  conformally  into  the  upper  half  of  the  t-plane  [8],  Figure  9. 


t 

M  /  +  M     (56) 


where  M  and  M.  are  complex  constants,  B,  C,  D  and  E  are  the  interior 
angles  (in  radians)  of  the  embankment  in  the  z-plane,  and  -1,  -P,  p 
and  1  are  points  on  the  real  axis  of  the  t-plane  corresponding  to  the 
respective  vertices,  B,  C,  D  and  E.  The  complex  constant  M.  corresponds 
to  the  point  on  the  perimeter  of  the  embankment  that  has  its  image  at 
t  =  0.  Substituting  the  values  of  B,  C,  D  and  E  into  equations  (56) 
leads  to: 
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B 


(a)       Z-  Plane 


} 


COD 


-P 


p      e 


(b)        t- Plane 
Figure  9.      The  Geometrical   Transformation 
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z  =  M  J  (  '\  '   xg  1  N  dA  +  ^  (57) 


where  N  e  —  ,  a  is  the  positive  single  of  the  slope  of  the  embankment. 
The  constants  M,  p,  and  M-,  are  determined  by  considering  the  corre- 
spondence between  the  z-plane  and  the  t-plane  at  points  D,  E  and  at  the 
center  of  symmetry. 
(a)  -  at  point  D,     t  =  p,    z  =  L  -  iH. 


K/     2        1 

Letting  /  I   g  — ^  )  N  dA  =  F   and  substituting  this  into  equation 

(57),   gives: 


L-iHoMF+J^  (58) 

(b)  -  at  point  E,  t=l,  z  =  L  +  H  cota. 

1 

Letting  /  !  -^ 2  )  N  dA  =  F  +  iF   and  substituting  into  equation 

0 
(57)  gives: 

L  +  H  cota  =  M(F1  +  iF2)  +  M^  (59) 

(c)  -  at  the  center  of  symmetry,     t  =  0,    z  =  -iH  and 

Then  ^  =  -iH  (60) 
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From  equations  (58)  and  (60)  one  obtains: 

M  -  §  (61) 

H^  =  -iH  (62) 

Substituting  these  values  into  equation  (59)  one  obtains: 

F      L 
F„  *  K 


or 


(63) 


F 

_2  H 

F   =  L 


Equations  (63)  represent  implicit  expressions  for  F,  F_  and  p  in  terms 

-  L     H 
of  s  or  z   . 

Figures  10,  11,  12,  13  and  Ik   show  the  solution  of  these  equations  for 

various  values  of  L/H  and  H/L.  Table  1  gives  the  values  of  F  and  p 

•Which  are  used  in  the  subsequent  illustrations. 

Substituting  equations  (6l)  and  (62)  into  equations  (57),  gives: 

t   2      1 

0 
This  equation  is  the  transformation  formula  which  maps  the  z-plane  into 

the  t-plane.  The  solution  to  equations  (63)  and  (6k)   are  evaluated  nu- 
merically on  the  IBM  709^  digital  computer  by  Simpson's  rule.  The  steps 
of  the  numerical  integration  are  obtained  by  trial,  so  that  the  numerical 
results  are  accurate  to  0.01  per  cent.  The  computer  programs  for  .the 
solutions  of  equations  (63)  and  (6k)   are  included  in  appendices  B  and  C. 

Operations  Connected  with  Conformal  Mapping  into  the  ilalf-Space 
Denote  the  region  under  consideration  by  S,  its  boundary  by  L,  and 
the  transformation  formula,  equation  (6k),   by  f(t),  then: 
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Integral    For  L/H  =  0.0l 
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|  -  f(t)  =  f(|  +   n)  (65) 

As  mentioned  above,  this  function  maps  S  into  the  lower  half  of  the 
t-plane,  i.e.,  into  the  half -plane  t)  >  0,  so  that  finite  points  corre- 
spond to  finite  points,  Figure  9.  Therefore,  straight  lines 
i)   =  constant,  lying  in  this  half -plane  correspond  in  the  z-plane  to 
some  lines  which  go  to  infinity  at  both  ends;  these  lines  will  be  de- 
noted by  I  . 

Similarly,  the  semi-infinite  straight  lines  |  =  constant  in  the 
lower  half  of  the  t-plane  correspond  in  the  z-plane  to  lines  x\   which 
begin  on  the  boundary  and  go  to  infinity  as  shown  in  Figures  9  (a)  an(i 
9  (b). 

Since  in  the  region  S  a  completely  definite  point  z  =  f  (|  +  ii\) 
of  the  z-plane  corresponds  to  every  point  (|,  rj)  for  r\  >  0,  the  quan- 
tities i   and  t)  may  be  conceived  as  curvilinear  coordinates  in  the  z- 
plane.  Due  to  the  nature  of  the  Schwarz-Christoffel  transformation, 
the  lines  £  and  r\   form  an  orthogonal  net  of  coordinate  lines.  If  z  is 
some  point  of  the  region  3,  Figure  15,  the  tangents  to  the  lines  £,  r\ 
at  z,  in  the  directions  of  increasing  |  and  tj,  will  represent  the  axes 
of  the  curvilinear  coordinates  at  the  point  z.  These  tangents  are  also 
denoted  by  |,  r,. 

Let  A  be  some  vector  starting  from  the  point  z  «=  f(£  +  iTj)  and  let 

A  ,  A  be  its  projections  on  the  x  and  y  axes,  respectively,  and  A  , 

x   y  % 

A   ,  its  projections  on  £  and  r\   axes,  respectively.  Then: 

A.  +  iA  =  (A  +  iA  )  e"i9  (66) 

£     T)   v  x     y'  v  ' 


k7 


where  G  is  the  angle  between  the  ox  and  o£  axes,  measured  from  the 
latter  anti-clockwise. 

If  the  point  z,  Figure  15,  be  given  a  displacement  dz  in  the  tan- 
gential direction  g,  the  corresponding  point  t  will  undergo  a  displace- 
ment d|  >  0  in  the  direction  g  of  the  t-plane.  Obviously, 

dz  =  eiejdz|  (67) 


whence 


19         dz  f'(t)  dg 

e    =fdzT=  iWtJikii 


and 


ie       f(t)  -ie        f(t)  ,,m 

e     =JFTtT]  .  e       ■  |f(t)l  (68) 


e2ie  m  _XM-  (69) 


Combining  equations  (66)  and  (69): 


Ag  ♦  iAs  -  (Ax  ♦  iAy)  1^*|-  (70) 


The  relationships  between  the  stress  in  the  oxy  and  o|tj  directions  are 
[35]: 


o\  +  a    =  a     +  a 

I  T    x    y 


a     -  0+   2iT„   =  (a  -  a  +  2ii     )  e 
tj    g      gij   v  y    x     xyy 


(71) 
2i0 


U8 


77  =  Constant 


0 


Figure  15.      The  Transformation  of  a   Vector   From 
The    Z-  Plane  into  The    t- Plane 
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Hence  the  following  relationships  hold; 


a    +  a 
x          y 

2 

-M 

a     - 

L  y 

0 

x 

+  2i 

O+O 
x+     y 

2 

»*4 

a     - 

.  y 

•  a 

X 

+  2i 

.] .-} 


e 
xy. 


It 
xyj 


2ie 


(72) 


■a    -  a 


ftP  * '  v]  *219} 


and 


a     +  a. 


i-th-v^J'-"'} 


-^♦*4[ 


a     -  o„   +  2iT  , 

n       5  ns 


]  *-2ie 


xy 


■  4 11^ + 1*,  J  »"; 


(73) 


The  Complex  Representation  of  Stresses  in  the  t-plane 
The  functions  0(z),  $(z),  4>'(z)  and  $'(z)  can  be  written  in  terms 
of  t,  in  the  following  forms: 


*(z)  =  Jf(t) 
$(z)  <=f[f(t) 


-•(t) 
=  f(t) 


(continued) 
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rw  -  m 


(74) 


Substituting  these  expressions  into  equations  (k2)   and  (^3),  gives: 
ov+  o^  =   2[0(t)  +  Ht)]  +YT7-   Im[f(t)] 

-U  Ee[*(t)]  +^7.  Im[f(t)]  (75) 


a  -  a  + 


2iV  =  j[f<t)  £^}  +  $(t)]  ♦  j=9S  7-  lm[f(t)] 


Combining  equations  (75)  and  equations  (71)  leads  to: 


o    +  a^^k  Re[*(t)]  +  ~I7Im[f(t)] 


a     ~  a,   +  2iT       ■  2 

n         i  xy 


fTtT         FTtJ 


1  ^  FTtT 


(76) 


Equations  (76)  require  determination  of  <J(t)  and  $r(t)  in  order  to  eval- 
uate the  stress. 

The  Complex  Representation  of  Displacements  in  the  t-plane 

The  displacement  in  the  z-plane  and  in  terms  of  the  complex  vari- 
able z,  is  given  by  equation  (5^).  Replacing  z  by  f(t)  in  this  equation: 
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2G(u  +  iv)  m   (3-J+n)*(t)  - 


m. 


TTtJ 


♦  *(t)  -  i|r(t) 


1-2M 


+  1^7  lm[f(t)]' 


(77) 


Applying  equation  (70)  to  equation  (77) »  the  components  of  the  dis- 
placement in  the  t-plane  can  be  written: 

f(t) 


2G(u  +  cd) 


(3-»*n)*(t) 


fTET 


TtT  -  t(t) 


i  =—«=  7  Im[f(t)]' 


f'(t 


(78) 


From  equations  (77)  and  (78),  the  following  expressions  hold; 


ucIg    R{ 


(3-4n)*(t)  - 


JM 


fTtT 


♦•(t)  -  +(t) 


+  i^7  lm[f(t)]2 


-w  4< 


=  ~     inl  (3-Uu)*(t) 


f(t) 


FTtJ 


♦•(t)  -  v(t) 


+  1^7  Mf(t)]2 


•-s  n 


(3-4u)«(t) ^L_  «.(t)    -  ir(t) 

FTtT 


^r^t)f]^l 


(continued) 
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\-    Im{  f  (3-M*(t)  -       ffi      ♦  •(t)  -  t(t) 

2G  I  L  fTTTT 


+  1^7  lm[f(t)]2 


\f$v^ 


Formulas  and  Propositions  for  the  t-plane 
As  mentioned  above,  the  expression  z  =  f(t)  transforms  the  region 

of  the  embankment  into  the  lower  half  of  the  t-plane,  which  is  bounded 

by  the  o\   axis  and  consists  of  the  points  t]  >  0. 

Assuming  that  the  stress  components,  which  satisfy  the  conditions 

of  continuity,  equilibrium  and  differentiability,  tend  to  zero  as  t  -*  », 

let  the  $(t)  and  $(t)  functions  be  represented  for  large  (t)  by: 


t    jt1+€ 


(30) 


t      |t1+ei 

where  c  and  c,  are  constants,  v  depends  only  on  t  and  v  -*  0  as  |t[  -*  <», 
and  e  is  a  positive  constant  which  is  arbitrarily  small.   Integrating 
the  above  equation: 

<>(t)  =  c  log  (t)  +  v  +  constant 

(81) 
\{r(t)  =  c,  log  (t)  +  v  +  constant 

in  these  expressions  one  definite  branch  of  the  multi-valued  function 
log  t  must  be  selected,  e.g.,  Log  t  +  i9,  where  0  (argument  of  t)  varies 
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from  -jt  to  jr.  In  addition,  the  functions  0(t)  and  ty(t)  will  be  har- 
monic in  every  finite  region,  contained  in  the  lower  half  plane  r\   >  0. 

Finally,  introduce  the  following  condition:  the  resultant  vector 
of  the  external  forces,  applied  to  a  segment  mn  of  the  o£  axis,  tends 
to  a  definite  limit  as  m  and  n,  Figure  16,  move  to  infinity  (m  towards 
the  left  and  n  toward  the  right).  This  condition  will  always  be  satis- 
fied  if  only  a  finite  part  of  the  boundary  is  loaded. 

This  last  condition  mathematically  means,  if  X. ,  X  are  the  com- 

i     n 

ponents  of  the  resultant  of  the  external  forces,  applied  to  mn,  then 
one  has  by  (32)  and  (38)  the  following: 


r  +  a  = .  ifafoi  7)  +  i  M*u±] 

5  n  l      c>x  5y       Ja 


(02) 


=  -  i 


♦(t)   +  f(t)  -|» 


The  negative  sign  has  been  chosen  on  the  right-hand  side  of  equation 
(82)  because  the  region  3  lies  on  the  right,  and  not  on  the  left,  as 
one  moves  from  m  to  n. 

When  calculating  the  increase  of  the  left-hand  side  of  equation 
(82)  for  the  transition  from  m  to  n,  one  may,  instead  of  moving  along 
the  straight  line  mn,  follow  any  curve  which  lies  in  S,  such  as  the 
curve  men,  Figure  16. 

If  m  and  n  lie  sufficiently  far  away  and  on  different  sides  of  o, 
then  equations  (80),  (8l)  and  (82)  give: 

X     +  iY     =  -il     c  log  r~-  +  cjri  +  c1  log  g-  -  c^iri  +  v  j      (83) 
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where  R  and  PL.  are  the  distances  of  m  and  n  from  0  and  v  is  arbitrarily- 
small  and  tends  to  zero  as  E  and  FL  tend  to  infinity.  In  order  that 
the  preceding  expression  will  remain  finite  for  any  arbitrarily  large 
R  and  R,  which  are  independent  of  each  other,  it  is  necessary  and  suf- 
ficient that : 

c  +  c^  ■  0  (8U) 


Under  this  condition  the  vector  (X  ,  X  )  of  the  external  forces,  applied 
to  the  whole  of  the  o|  axis,  will  be  given: 


X^  +  iY^  =  -i[*(t)  +  f (t)   ^|   7T0j_=  «(c  -  cp   (85) 
It  follows  from  (81+)  and  (85)  that: 


Xe  +  iY  X  -  iY 


Hence,  for  large  |t|: 


Xe  +  iY 

#(t)  =  -  Act-3  + 


.1+6 


«*'  "  ^53T*  +  ^  <8?> 


Xt  -  iY      u2, 
2nt2.     |t2<1+€' 
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Xe  +  iY 
<t»(t)  =  -  »  %   "  log  t  +  v  +  constant 


X   -  iY 
ty(t)  =  ■  ■  i  3  log  t  +  v  +  constant 

The  quantities  Y  =  N(t)  and  X  =  T(t)  have  to  be  found  on  the 
axis  of  ^  as  a  function  of  the  abscissa  £ .  On  the  basis  of  equation 
(87),  for  large  t: 


,.<i),     T  =  <|) 


(88) 


where ; 


l<J><W 


These  conditions  are  necessary  conditions  in  order  to  apply  the 
Cauchy  integral  formula  [25],  which  is  introduced  below. 
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BOUNDARY  CONDITIONS 

Boundary  Conditions  in  the  z-plane 
From  equation  (75) 


a     +  a  a     -  a    +  2iT 

N  +  iT  =   _2i_y.   +  JL ^ SL 


or 


N  +  iT  =  3>(t)  +  *(t)  +  z$'(z)  +  iy(z)  +  ?Im[f(t)]       (89) 


Denoting  the  boundary  of  the  embankment  by  L,  equation  (89)  becomes: 


;  N  +  iTj  ■     *(z)   +  $(z)  +  z*'(z)   +  f(z)   +  7lm[f(t)]L      (90) 


or 


ll  +  iTjL=  12  Re[o(z)]   +  ztf'(z)   +  y(z)   +  y  Im[f(t)]  jL     (91) 


where : 


z  =  x 

z  «s  x  -  i 

z  =  x  -  iH 

z  =  x  -  i 

z  =  x 


L  +  H  cotq  +  x 
cota 


L  +  H  cota  -  x 


cota 


for  x  <  -  (L  +  H  cota) 
for  -(L  +  H  cota)  <  x  <  -L 
for  -L  <  x  <  L 
for  L  <  x  <  L  +  H  cot 
for  L  +  H  cota  <  x 


Boundary  Conditions  in  the  t-plane 
Ft  om  equat  i  ons  ( 76 ) : 
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a    +  a„       a    +  a„   +   2i-r, 
2    -I-       2 


N  +  iT«  3_1+_ji — I La 


N  +  iT  «  *(t)  +  «(t)  + 


"?(ty   »'W    +     *'W  f(t)l 

L     FTtJ      FTtT    J 


f'(t) 


(92) 


The  boundary  of  the  embankment  in  the  z-plane  corresponds  to  the  axis 
t)  =  0  in  the  t-plane.  Substituting  this  in  equation  (92)  gives: 


n  +  iT  »  *(i)  +  *(D  +  ff(i)    ~jjj     +     f'^e)    $(0 

L  7TTT  ?Ti  1  J 


fTTT         fTFT 


+  Ti-|la[f(6)]  +  ^r|    lm[f(6)]       f'^}  (93) 


fTIT 


Because  the  boundary  is  free  of  load,  N  +  iT  =  0,  and  equation  (93) 
can  be  rearranged  to  yield: 


*(l)  +  *(|)  + 


--L.Zlm[f(e)]  -i=^|lm[f(e)]e 


2ia 


(9»0 
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where : 

a 


for 

-     CO     <     | 

<  -1 

for 

- 1  <  1 

<  -0 

for 

-  P<£ 

<    3 

for 

p  <  1 

<  l 

for 

1  <  I 

<    CO 

and 


Im[f(|)] 


0         for 

--eo  <  | 

< 

-i 

■Hp0   for 

-  1<  g 

< 

-p 

-H         for 

-  P  <  1 

< 

p 

=5^ir» 

P     <  1 

< 

i 

0         for 

i    <  i 

< 

CO 

Thus,  the  boundary  conditions  are  as  shown  in  Figure  17.  These 
conditions  will  be  called  "Fictitious  boundary  stresses",  because  they 
arise  from  the  weight  of  the  embankment. 

It  is  convenient  to  separate  the  real  and  imaginary  parts  of  the 
"fictitious  boundary  stresses",  i.e., 


and 


TwE  -  Imj-  lEr  I  ^f(S)]  -  fer  i  im  tf(|)l  e2iG}    (96) 
where  BL_  and  T,„  are  shown  in  Figure  17  (b)  and  Figure  17  (c), 
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•00— A  By^ 


F-co 


(a)     The    Boundary  of  the  Investigated    Problem 


(b)     The    Real   Part    of   "Fictitious  Stresses"    on 
the    Boundary    in    the    t  -Plane,  Nwe 


r^ 


~c         D- 


oo  —  -I  -p  p     £       I  -—  oo 


(c)     The    Imaginary   Part  of  "Fictitious   Stresses" 
on  the   Boundary    in   the    t-  Plane,  TWE 


Figure  17.        "Fictitious  Boundary    Stress"    Due    to 
the    Weight   of   the    Embankment 
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respectively. 

Incorporating  these  notations  into  equation  (9^)»  leads  to 


l\lE  +  iTWE  =  *<*>   +  «(6)   + 


fHTT         fTTT 


(97) 


and  by  taking  i  to  be  -  i 


VB-^-tCe) +#(*)+   f(i)4^f +-f[||-$(o_ 


These  equations  lead  to  the  solution  of  the  problem,  as  shown  below. 
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THE  SOLUTION  OF  THE  PROBLEM 


Introduction 


In  the  following  section  the  use  will  be  made  of  so-called  Cauchy 
integrals.  A  systematic  study  of  the  properties  of  these  integrals  may 
be  found  in  reference  [39]»  hut  for  convenience,  the  essentials  will  be 
given  in  the  next  section. 

Cauchy  Integral 

Let  h(£)  c  h..(£)  +  h?(|)  may  be  some  complex  function  given  on  the 
boundary  L  of  the  embankment.  Assume  that  h(|)  is  a  finite  and  inte- 
grable  function  in  an  ordinary  Riemann  sense. 

The  integral  of  the  form 


taken  over  L  with  |  some  point  on  the  boundary  L,  is  called  a  Cauchy 
integral  [39] •  It  is  assumed  that  the  point  t  does  not  lie  on  the 
boundary  L. 

The  above  integral  represents  a  function  of  the  complex  variable 
t  throughout  the  entire  plane  with  the  exception  of  the  points  L.  De- 
note this  function  by  H(t),so  that: 
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The  function  H(t)  is  harmonic  inside  the  region,  and  if  h(t)-»  €  as 
t  -►  oo,  then  H(t)  -*  e  [39]>  where  e  is  some  finite  value. 

From  the  properties  of  Cauchy  integrals,  the  following  theorems 
can  be  stated  [39]« 

1  -  If  the  function  H(t)  is  harmonic  and  continuous  inside  the 
region  and  H(t)  I.     ■  0.  Then 

2«i  /  E[  ft   "  H^      inside  the  region  (100) 

and 

■—  I   H^d|  ■  0        outside  the  region  (101) 

2  -  If  the  function  H(t)  is  harmonic,  continuous  outside  the  region  and 
H(t)|.     =  0.  Then 

1 1   =  CO 

2nl  /  6  -  t  ~   °        inside  the  region  (102) 


and 


~r  /  Hffid|  =   -  H(t)  outside  the  region  (103) 


l*- 


2«iJ   I  -  t 


The  expressions  (100),  (101),  (102)  and  (103)  are  called  the  Cauchy 
formulas . 

Determination  of  the  Stresses  by  Using  Cauchy  Integral 
It  was  shown  above  that  the  mapping  formula  is  given  by: 


*,   2 
z   L    I  A  -  1 


H   HFJ  L  • 


N  dA  -  i  =  f(t)  (lOU ) 


0    A  "P 


6U 


When  the  point  t  moves  along  the  |  axis  from  left  to  right,  the  corre- 
sponding point  z/H  moves  along  the  boundary  likewise  from  left  to  right. 

It  is  readily  verified  that  the  equation  (10U)  maps  the  region 
under  consideration,  into  the  half -plane  ij  >  0. 

Repeating  equation  (97)  for  reference, 

*(i)  +  ^TiT  +  ffTiT      *'(§)   +  -£-f*f-  g(e)1  -  n..„  +  t      (105) 

In  this  equation,  the  unknown  functions  <E>(t),  f(t)  which  are  harmonic 
in  the  lower  half  of  the  t-plane,  satisfy,  on  the  basis  of  (87),  the 
conditions : 


•<*>(% -.-  w 


••Wit---    rh  (lo6) 


*<*>!*...-  iTf 


Therefore,  on  the  basis  of  equations  (100),  (101),  (102)  and  (103) 
the  Cauchy  integral  can  be  applied  on  equation  (105). 

The  function  f(|~) — : — *** —  determined  by  (105)  is  the  boundary  value 
of  the  function  g(t)  ■  ■  ■*■  ■'■   harmonic  in  the  lower  half  of  the 

frTty 

t-plane  and  vanishing  at  infinity,  and  the  function  *Sj  f(|)  is 

f  •  (t) 

the  boundary  value  of  the  function  — ■  ■  v  '   $(t),  harmonic  in  the 

FTty 


65 


lower  half  of  the  t-plane  and  vanishing  at  infinity.  Applying  the 
Cauchy  integral  formula  to  equation  (105),  after  taking  the  complex 
conjugate  form  of  it,  gives: 


"TCO  "TOO  "TOO 


2rciJ      |-t       2niJ      g-t       2rtiJIU;frtT)|- 


2*iJ   f'(£)  I  -  t  Q§   2*ij     |-t  UU° 

—OO  —CO 

% 

where  t  is  a  point  of  the  lower  half  of  the  t  plane;  noting  that  <5(fj) 
is  the  boundary  value  of  <*>(t),  harmonic  in  the  lower  half  of  the  t-plane 
and  vanishing  at  infinity,  and  that  <I>(|),  f(|)  f , K(  and  ftrfH  f(0> 

are  the  boundary  values  of  0(t),  f(t)  y| ■)■■■(  and  „t  ).i   $(t),  harmonic 
in  the  upper  half  of  the  t-plane  and  vanishing  at  infinity.  From  the 
properties  of  the  Cauchy  integral,  it  is  clear  that  the  first,  third 
and  the  fourth  integrals  vanish  and  the  second  integral  is  equal  to  $(t), 
Therefore : 

0(t)  =  SSj    g  -  t    d*  (108) 


The  function  $(t)  is  now  easily  determined  from  its  boundary  value  given 
by  equation  (105): 

f(t) .  .E(t)  m .  «t)  4^_+ ^_  &p**&« 

—00 

(109) 
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where 


«<*>  ■  wj 


LA2  -  p' 


dA  +  i 


(no) 


e2i0!   (111) 


^-^  =  -^1*^)3-  ^i^Wl)]--^  (a«) 


Substituting  these  values  into  equations  (108),  leads  to; 


-1 


...     1,  7g   ,    .     d& 

+  2rtiJ    2(l-n)        |    -  t 


+  53[ J  2(i-M(i-p)      fTT    d| 


-2sr~J  2(i-,o(i-p)    i  - 1  d^ 

-1 


+  £f^  /■    yH  .   .  JiL_ 

2rti        /    2(1-|.0        |   -  t 


-2ia3 


'  ^ssr1  J  2(i-Jhi-p)  *  K"t  d* 


(113) 


where : 
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°L 


a2  =  o 


°3"  N 


Integrating  and  elaborating: 


*(t)  _       i 

7H  -  -  4n(l-n)(l-pJ 


log  |g  +  t  log  |g  ♦  1-b] 


i    ,   t-B 
U«(l-n)  log  t+p 


Llog  t=£  - t  log  H  - 


i  e-2iaL 
i4rt(l-u)(l-6) 


log  £f  +  t  log  |±f  +  1  -  6 


i  e-2ia2  ,   t-B 


i  e-2i0^ 
Uk(1-h)(1-6) 


log^  -  t  log  |=|-1+  b]     (11U) 


Differentiating  this  equation  gives 


tt'(t)         i      [     1-p      ,   t+S_     t(l-B) 

7H  =  "  J+jt(i-n)(i-ej  I  (t+i)(t+6)    g  t+i   (t+i)(t+e). 


2n(l-n) 


Lt 


2   ,,2 


(continued) 
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TwT 


i  1-6  .       t-1  t(i-p) 

1-H)(1-P)  L(t-l)H-6)   -  l0S  tip  -   (t-l)(tip)J 


i  e"2iaL     ("  1-p  ...        t+6  t(l-B) 

c(w)(i-p)  L':t+i)(t+p) +  log  t+i    it+i)(t^. 


i  e-2ia2     [        p       1 


i  e"2iQ3     T        1-P  t        t-1  t(l-6)   ]    /nncx 

«(1-H)(1-P)  L(t-l)U-6)  "  l0g  tip  "  (t-lj(t-p)J   (u-5) 


Combining  equations  (109),  (110),  (ill),  (llM  and  (115),  gives: 


3felm  ■  -£&(:  ,   ,    *        . 

7h  FTtT  I  un;i-n)(i^) 


i-P  -L     T  t+p 

(t+iKW)     log  t+i 


(t+i)(t+6)j   ■  2^r^ry 


2-2 
L   t*  -  p"  J 


4«(i-i.t)(i-e) 


1-P        .     t-i       td-e) 

(t-l)(t-p)  - los  tip  -  (t-l)(t-p). 


i  e"2iaL       T         1-p  .  .       t -W5  t(l-p) 

Wi-n)d-g)  L  (t+iiU+p) +  log  s$ +  (t+i)(t^). 


i  e"2ia2 


2it(l-n) 


■t&   -   p*"  -• 


i  e"2ia3     [  1-P  !        t-1  t(l-p) 

"  Mi-n)U-p)L(t-iMt-p)  -  iog  tip  -  (t-ij(t-p). 


1 


m 


I  -   lht(l-n)(l-p) 


iog£34  *  iogUI+l' 

(continued) 


Urt(i-n)  XUB  t+p 


log 


t-£ 


Mi-t-OU-P) 


t-i  .  .   t-i 

log  5-  -  t  log  ^  -  1  + 


i  e-2iaL 

Mi-u)(i-p) 


lag  gf  +  t  log  gf  +  i 


i  e-2ia2  ,    t-B 

s^D — los  t^ 


.  -2ia_ 

i  e    3 

4«(l-n)(l-p) 


^H-tlogU-1  +  t}} 


,  f;(t)j        i       r 
ml*  ^(l-njd-p)  L 


log  ^  +  t  log  — E  +  1 


U(i-n)    logl^ 


l!rt(l-H)(l-f) 

i  e+2ial 


log  |=i  -  t  log  |=i  -  1  +  p 


i  e  "*"!   .   t+e  .  .  .   t+e  .  . 
-  Mi-uKi-p)  Llos  t3  +  t  log  ^  +  i  -  p_ 


i  e+-ia2  ,   t-g 

U(l-;0    l0g  t+1 
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.  +2ia, 

i  e    3 

4n(l-nKl-P) 


,   t-1   ,  .   t-1   .  . 


(116) 


Combining  equations  (lln),  (115)  and  (116),  gives: 
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-Zicu 


5;,2>L.^^ 


?H 


L"  ti*(i-n)(i-B) 


loe  IS  +  *  lo*  M  +  1  -  p] 


ar(i-n)    logl+f 
Ki4e-2ie3) 


-      t-1      .  -      t-1     -   . 
l0s^-tl0gt^-1  + 


f(t)-g(t) 

f'(t)  J 


+  Re! 
L 

+  (t+l)(t+f?)/  TTJIUT  ^        2  2  y 


Mi-u)(i-e)    V(t+i)(t+p)    +  xo&  t+1 


£J±1 


i(l+e"2ia3)  ,  1-B  .       t-1  t(l-p)    {] 

-  U«(i-,x)(l-p)     '^(t-l)(t-3)   "  l0S  tip  "   (t.lj(t.p);J 


i(l+e       °l)  /'         t+6   .    ,    ,        t+S    ,    . 
hn(U)(l^)     [Kl0&  t+1  +  *  log  t+T  +  X  " 


i  ,        t-6 

2Trr^T  1o^ 


i(l+e  3) 


t-1 


Mi--?(i-P)Alog  S  ■ t  log  t^  - x  +  V 


+  "TTt7 


i(l+e+2i°l)  /'  t+6  A   .    .        t+B       . 

"  W(l-^)(l-pj  Vl0S  tS  +  *  l0g  t+1  +  X  -  py 


i  ,       t-S 

2n(l-M)       g  t+p 


i(l+e+2iG3)  (.        t-1 
Ml-nKl-P)  'vl0g  t^  - 


(continued) 
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+  Im[f(t)] 


(117) 


ax  ,Eej  ,  i(^;2iap 
^  - 2  Re  r  i»x(U)(i-p) 


log  |i|  +  t  iog  |if  +  i 


i  .        t-B 

2«(l-n)   l0g  t$ 


i(l+e         j 
U«(l-u)(l-p) 


,        t-1       .    ,        t-1       ,    . 

log  tip  -  t  log  tip  "  L  + 


Re 


f(t)-g(t)T  i(l+e-2l0i;   ,  1-p  .    .        t+6 

Kt)  I"  Mi-p)(i-p)    ^(t+iKt+p)  +  los  ua 
(t+i)(t+p)  /  7{TTiIT  V    t2  -  p2  ^ 

i(l+e'2ia3)  ,  '         1-p  .       t-1  t(l-p)    \\ 

-  j4«(i-,o(i-p)'  ^(t-i)(t-p)  - log  tip  -  (t-i)(t-P);J 


f'(t)     r  i(l+e""ial)/(       t+e^..       t+B  ^  n        _\ 

i  ,        t-| 

-  2«(l-n)  l0St^ 

i(l+e~2ia3)  -  ,       t-1       .    _       t-1       .    .    _  \\ 

-  uli-Tod-p)  vlog  tip  -  * log  rip  - x  +  M 

f*(t)     [  i(l+e+2l<XL     ■  t+p  ^   .    ,        t+p    .    . 

?TO~  L  "MMN'v log  t3  +  *  log  t+f  +  L  -  V 

"  2«(l-u)    l0g  fel 


(continued) 
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i(l+e2ia3) 


t-1   ,  .   t-1 

logH-tlog^-1  + 


+  TB-  Im[f(t)] 


(118) 


_xy 
7H 


=  Im. 


f(t)  -  g(t) 
f'(t 


i(l^"2ial)  (        1-fi 
Un(l-n)(l-p)  ^(t+l)(t+0) 


,    t+fi  .    t(l-6)  \    i     /   6 
log  1+1  +  (t+l)(XfO>  ^T=TTT  ^2 


i(l+e-2i03)  /     1-p         t-1      t(l-p) 
Un(l-n)(l-p)  ',  (t-l)(t-p)  -  i06  t-p  "  ( t-1 j (t-p) 


f'(t)  [      i(l+e"2iai)  /'     t+6  .  .  ,    t+p  ,  . 

-  tW     '  Mi-»0(i-p)     v log  t+f  +  t  log  t+f  +  x 


2fl(i-fo  logfe§ 


i(H-e"'1'ia5)  ,■   .        t-1        .    .        t-1       .    .    .    \ 

itrt(i-;)(i.p)/'vi°g  tip  -  t ios  t^  - x  +  p  j . 


m[-M^(Si     (^Ig^los^.X-p 


i         ,       t-p 
2«(1-M)        g  t+p 


i(l+e+2ia3)     n       t-1       .   .       t-1       ,    .  QN 

14     (1-"(1-P)       I  lQS   ^-tl0^-1    +    P 


-(119) 


Determination  of  the  Displacements 
The  components  of  the  displacements  are  given  by  equation  (77), 
repeated  here  for  reference: 
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SUu  +  v]  =   (3-%)*(t)   -  $Q  *»(t)   -  *(t)   +  i  ifSii  lni[f(t)]2   (120) 


Substituting  the  integral  of  equations  (llU)  and  (ll6)  and  the  complex 
conjugate  of  equation  (llU),  into  equation  (120),  leads  to: 


^   +  ^  =  "  i(3"tnSl^Wl-u)"2iai][  (l-P)(t+P)log(t-*3) 


^(i+h) 


2  2 

+  (t  +  P)"  log  (t-ff)  -  -&-±H    log  (t  +  1) 


(  1  j;  F  )  (l  +  2t  +  p)  -  p  log  p  +  §-  ] 


£-i] 


i(3-V)[l+(l-2n)e"3lQ33 
Ml-P)U-u) 


(p-l)(t-p)log  (t-p) 


2  ^ 

+  ^  log  (t-p)  -  &giL  log  (t-1) 


+  -fcil  (2t-l-p)  -  p  ( |  +  l)  log  p  +  I  p 


V2   V 


+  in(l-p)  to)  +  J 


-  i(gf^  [(t-p)  log  (t-p)  -  (t+p)  log  (t-ffj) 

] 


+  2  p  log  p  +  inp 


(continued) 


7U 


-  ^fot&T11  [(i-e)(*+p)  iog  ™ 


+  ii|&L  log  (t+p)  -   (t+1)fa  log  (t+D 


+  (  ij&  )(l+2t+p)  -  p  log  p  +  |-  log  p  +  g-  -  | 


i[l-Kl-2|i)e"2ia31 
"     im(l-p)(l-n) 


(p-l)(t-p)  log  (t-p) 
+     ^fr     log  (t-p)  -  fi£L  log,(t-l) 


+  Ifijjii.  (2t-l-p)   -  p  (|  +  l)  log  p  +  I  p< 
+  i«(l-p)  (2±§J+  i 


b[^ 


p)  log  (t-p)  -  (t+p)  log  (t+p) 


+  2p  log  p  +  irtP 


f(tl»(t)(tT(T)  +  i  i^T  m[f(t)]2      (121) 


The  horizontal  displacement,  u,  is  the  real  value  of  equation  (121) 
and  the  vertical  displacement,  v,  is  the  imaginary  value  of  equation 
(121). 
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RESULTS  AND  DISCUSSION 


Introduction 


The  stress  distribution  within  and  under  long  linear  elastic,  iso- 
tropic homogeneous  embankments,  continuous  with  the  underlying  material 
of  similar  properties,  has  been  investigated  for  the  cases  shown  in 
Table  1  for  Poisson's  ratio,  p.  =  0.3  and  0.5.  These  results  were  com- 
pared with  those  obtained  by  the  normal  loading  approximation,  viz. 
Figure  6  (b).  All  the  results  below  show  only  the  effect  of  the  embank- 
ment weight.  Thus,  at  depths  below  the  base  of  the  embankment  (y/H  =  0) 
the  material  is  assumed  weightless.  The  effect  of  the  medium  weight 
can  be  superimposed  upon  these  values  to  give  the  total  stress  acting 
at  a  point. 

Analysis  of  the  Vertical  Normal  Stress  Distribution 
Figure  18  shows  a  typical  result  of  the  vertical  normal  stress 
contours  in  dimensionless  form,  a  JyYL,   for  an  embankment  with  a  =  ^5°, 
L/H  =  3,  and  Poisson's  ratio,  \i   =  0.3.  The  dashed  lines  in  this  figure 
show  stress  contours,  in  terms  of  a   /7H,  for  the  usual  normal  loading 

y 

approximation  corresponding  to  this  embankment.  These  indicate  that 
the  vertical  normal  stresses  produced  in  the  foundation  material  below 
the  elastic  embankment  are  generally  smaller  than  those  computed  for 
the  normal  loading  approximation. 

The  stress  distribution  due  to  the  normal  loading  approximation 
*  Experiments  suggest  \x   =  0.3  is  a  reasonable  value  for  unsaturated  soils, 
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Figure  18.     Contours  for  Vertical   Stress 
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is  independent  of  Poisson's  ratio;  the  stresses  due  to  the  elastic  em- 
bankment are  dependent  upon  u.  However,  the  vertical  stresses  are  rela- 
tively insensitive  to  its  magnitude;  changing  \x   from  0.3  to  0.5  changes 
the  vertical  stress  at  a  point  by  less  than  five  per  cent. 

At  the  base  of  the  embankment  (y/II  «*  0)  there  is  a  sharp  change  in 
slope  of  the  contour  lines.  This  results  from  the  difference  in  the 
way  that  a  weightless  body  is  stressed  under  an  elastic  body  force  ap- 
plied to  its  surface.  The  expression  for  the  vertical  normal  stress, 
equation  (117)  has  two  parts.  The  second  part  of  this  equation  is  a 
linear  function  of  the  unit  body  weight  and  position.  It  is  expressed 
by  7lm[f(t)].  In  the  foundation  y   is  zero.  Therefore,  as  the  stress 
contour  passes  over  the  interface  line  a  linear  term  is  dropped  from 
the  equation  for  stress.  This  results  in  the  observed  sharp  change  in 
the  contour  lines. 

The  effect  of  L/H  ratio  on  the  vertical  stress  along  vertical  sec- 
tions through  the  center  line  of  the  embankment  and  the  toe  of  the 
slope,  is  illustrated  in  Figure  19,  for  a  =  1+5°,  n  ■  0.3.  The  figure 
is  a  composite  diagram  showing  the  embankment  schematically,  and  the 
magnitude  of  the  vertical  stress  at  each  section  as  a  function  of  depth. 
Clearly,  the  L/H  ratio  has  a  pronounced  effect  on  the  distribution  of 
vertical  stress.  As  L/H  decreases,  the  stress  decreases.  Furthermore, 
a  smaller  L/H  ratio  produces  a  more  rapid  dissipation  of  stress  with 
depth,  and  elevates  the  horizontal  section  where  the  vertical  norma,] 
stress  becomes  practically  uniform. 

The  dashed  lines  show  the  vertical  normal  stress  for  the  normal 
loading  approximation  equivalent  in  shape  to  the  embankment  for  which 
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y/H 


(a)  At   Centerline  (b)  At  Toe  of  Slope 

Figure- 19.     Distribution    of      Vertical    Stress 
Along     Vertical    Sections     for 
Varying    L/H     Ratios 
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L/H  =1.  It  is  clear  that  the  vertical  stress  beneath  the  embankment 
for  the  corresponding  elastic  embankment  is  smaller. 

Figure  20  shows  the  distribution  of  vertical  normal  stress  along 
the  base  of  the  embankment  for  p.  =  0.3,  four  values  of  a,  and  several 
embankment  shapes  shown  schematically  in  the  figure.  The  curved  solid 
lines  represent  the  distribution  of  stress  against  the  base;  dashed 
lines  show  the  distribution  of  stress  assumed  in  the  usual  normal  load- 
ing approximation.  It  is  clear  that  the  stress  distribution  is  much 
more  uniform  under  the  elastic  embankment  than  is  ordinarily  assumed. 
The  difference  becomes  especially  apparent  as  the  L/H  ratio  of  the 
steeper  embankment  decreases.  Moreover,  the  magnitude  of  the  stress 
under  the  central  zone  of  the  elastic  embankment  is  less  than  that  shown 
by  the  dashed  curves.  Again  the  effect  is  enhanced  for  narrow,  steep 
embankments  (for  a  =  ^5°,  and  L/H  =  0,  the  vertical  normal  stress  is 
only  65  per  cent  of  that  usually  assumed.) 

Because  there  are  no  shearing  stresses  on  the  axis  of  symmetry  of 
the  embankment,  the  conditions  of  equilibrium  require  that  the  sum  of 
the  vertical  normal  stresses  acting  on  the  interface  must  be  equal  to 
the  weight  of  the  embankment.  Therefore,  the  area  under  corresponding 
dashed  and  solid  curves  must  be  the  same.  Hence  the  difference  betv/een 
these  curves  becomes  less  pronounced,  at  least  near  the  central  portion 
of  the  embankment,  as  the  L/H  ratio  increase.  However,  near  the  outer 
edge  of  the  embankment,  the  stresses  are  still  significantly  larger  on 
a  proportional  basis,  than  indicated  by  the  normal  loading  approximation. 
Thus,  for  embankments  with  moderate  L/H  ratios,  the  normal  loading  ap- 
proximation leads  to  larger  estimates  of  differential  settlement, 
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Figure  20.  Distribution      ot     Vertical     Normal     Stress,  a,  fyW     On      The    Base 
of     The     Embankment     For     Varying    oc.  And  L/H     Ratios 
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assuming  one-dimensional  compression,  than  would  be  computed  by  the  meth- 
od presented  herein. 

Analysis  of  the  Horizontal  Normal  Stress  Distribution 
The  horizontal  normal  stress  contours,  a  /yE,   are  shown  in  Figure 
21,  for  a  =  U50,  L/H  =  3  and  u  =  0.3  The  dashed  lines  are  contours  de- 
termined from  the  usual  normal  loading  approximation.  The  figure  shows 
that  the  maximum  horizontal  stress  occurs  within  the  body  of  the  embank- 
ment and  decreases  with  increasing  depth.  In  the  foundation  material  in 
the  vicinity  of  the  elastic  embankment,  o  /7H  is  less  than  half  of  that 
usually  assumed,  because,  the  normal  loading  approximation  neglects  the 
shear  stress,  t  /7H,  transmitted  from  the  embankment  to  the  foundation. 

The  effect  of  the  ratio,  L/H,  on  the  horizontal  stress  along  ver- 
tical sections  through  the  center  line  of  the  embankment  and  the  toe  of 
the  slope,  is  illustrated  in  Figure  22  for  a  =»  1+5  °  and  u  =  0.3.  As  the 
embankment  becomes  narrower  (L/H  <  l),  the  stress  is  actually  negative 
at  some  points  below  the  center  line.  That  is,  the  embankment  causes  a 
reduction  in  horizontal  stress  at  these  points.  However,  the  horizontal 
normal  stress  within  the  embankment  at  the  center  line  increases  as  L/H 

increases  and  becomes  larger  than  b*L  for  L/H  >  3.  This  is  due  to  the 

1-u        - 

shear  stress  acting  along  the  base  of  the  embankment. 

The  dashed  line  shows  the  stresses  determined  from  the  normal  load- 
ing approximation  for  L/H  =  1.  The  stress  is  larger  than  that  due  to 
the  elastic  embankment  at  all  depths.  In  fact,  in  the  vicinity  of  the 
embankment  it  is  more  than  five  times  as  large  under  the  center  line 
and  twice  as  large  under  the  toe.  In  contrast  to  the  elastic  embank- 
ment, the  normal  loading  approximation  does  not  produce  negative 
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a=45° 

L/H»3 
^  =  0.3 


Figure  21.  Contours      for       Horizontal        Stress 


83 


(a)  At  Centerline         (b)  At  Toe  of  Slope 

Figure 22.  Distribution     of    Horizontal     Stress 
Along     Vertical     Sections    for 
Varying    L/H     Ratios 
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horizontal  stress  at  any  depth.  The  reason  for  this  difference  becomes 
apparent  when  the  shear  stresses  transmitted  by  the  embankment  to  the 
foundation  material  are  considered.  This  is  discussed  below. 

The  effect  of  Poisson's  ratio  on  the  horizontal  stress  is  illus- 
trated in  Figure  23.  This  figure  shows  the  horizontal  stress  along  ver- 
tical sections  through  the  center  line  and  toe  of  the  embankment  for 
a  =  U50  and  p.  =  0.5.  The  dashed  line  shows  the  stress  due  to  the  nor- 
mal loading  approximation  for  L/H  =  1.  Comparison  with  Figure  22  in- 
dicates that  a  change  in  Poisson's  ratio  from  0.3  to  0.5  changes  the 
stress  at  shallow  depths  below  the  central  portion  of  the  embankment  by 
as  much  as  a  factor  of  two.  The  difference  decreases  as  the  L/H  ratio 
increases.  The  influence  of  u  is  less  pronounced  below  the  toe  than 
below  the  center  line. 

Analysis  of  Horizontal  and  Vertical  Shear  Stress 
Contours  of  horizontal  and  vertical  shear  stress,  1     /yE.   are  shown 
in  Figure  2k,   for  a  =  45°,  and  L/H  ■  3.  The  solid  contours  are  for 
Poisson's  ratio  of  0.3.  The  long  dashed  contours  are  for  \jl   b  0.5,  and 
the  short  dashed  lines  are  for  the  normal  loading  approximation.  The 
figure  indicates  the  existence  of  horizontal  shear  stresses  within  the 
body  of  the  embankment,  increasing  to  a  value  at  the  base  near  the  toe 
of  the  slope,  in  excess  of  O.27H.  However,  the  maximum  value  of  hori- 
zontal shear  stress  (approximately  O.37H)  occurs  below  the  base  of  the 
embankment . 

Like  the  horizontal  normal  stress,  the  shear  stress,  1  /yYl,  is 
affected  markedly  by  the  magnitude  of  Poisson's  ratio.  However,  the 
effect  observed  depends  uon  the  position  of  the  point  considered, 
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(a)  At  Centerline       (b)At  Toe  of   Slope 

Figure  23.  Distribution   of    Horizontal    Stress 
Along    Vertical    Sections   for 
Varying    L/H    Ratios 
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y/H 


12 


13 


Figure  24.    Contours    of      Shear    Stress,  ~C: 
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relative  to  the  base  of  the  embankment.  In  the  zone  below  the  embank- 
ment to  a  depth  of  y/H  equal  approximately  two  to  three,  the  shear 
stresses  in  the  incompressible  material  (u  =  0.5)  are  less  than  for  the 
case  in  which  u  =  0.3.  At  greater  depths,  the  reverse  is  true.  The 
shear  stress  determined  from  the  normal  loading  approximation  is  less 
than  that  for  either  u  above  a  depth  factor  of  approximately  three  to 
five,  and  more  at  greater  depths.  The  magnitude  of  this  effect  depends 
upon  the  horizontal  location  considered,  as  shown  in  the  figure. 

The  normal  loading  approximation  assumes  that  there  is  no  shear 
stress  at  the  base  of  the  embankment.  Figures  2k   and  25  indicate  that, 
for  the  elastic  embankment,  this  assumption  is  not  reasonable.  Figure 
25  shows  the  horizontal  shear  stress,  x   /yH,  at  the  base  of  the  embank- 
ment for  n  =  0.3,  four  values  of  a,  and  a  variety  of  embankment  shapes 
shown  schematically  in  the  figure.  The  horizontal  shear  stress  is  zero 
at  the  center  line,  as  required  by  symmetry,  and  reaches  a  maximum  near 
the  toe  of  the  slope.  The  magnitude  of  the  maximum  and  its  location 
depend  upon  a  and  the  embankment  shape.  As  L/H  decreases  for  a  given 
a,   the  maximum  t  /7H  increases,  and  moves  closer  to  the  toe  of  the 

^y 

slope.  The  magnitude  of  the  increase  is  slight  for  a  =  15  degrees,  but 
becomes  more  significant  as  a  increases.  Note  that  a  maximum  t  /7H 
in  excess  of  O.k   implies  that  the  horizontal  shear  stress  at  the  base 
of  a  forty  foot  high  embankment  may  be  greater  than  one  ton  per  square 
foot  (unless  the  shear  strength  of  the  material  is  such  that  failure  is 
induced) . 
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Figure  25.       Distribution     of      Shear      Stress,    T„/;rH      On      The     Base 
of     The     Embankment      For     Varying    a,  And   L/H     Ratios 
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Analysis  of  the  Maximum  Shear  Stress 
It  is  often  useful  to  consider  whether  the  maximum  (i.e.,  princi- 
pal) shear  stress,  t  /7H,  at  any  depth  beneath  the  embankment  exceeds 
the  available  shear  strength.  Thus,  it  is  desirable  to  know  the  magni- 
tude and  distribution  of  maximum  shear  stresses  due  to  the  embankment. 
Contours  of  1       /yE  are  shown  in  Figure  26  for  the  embankment  section 
of  Figures  18,  21  and  2U.  Note  that  the  magnitude  of  t    transmitted 
from  the  embankment  to  the  foundation  material  is  approximately  O.257H 
at  the  base  of  the  embankment  in  the  vicinity  of  the  toe.  However,  the 
largest  shear  stress,  O.337H,  occurs  beneath  the  center  line  at  y/H  =  1.? 
It  is  also  interesting  to  observe  that  within  the  embankment,  the  maxi- 
mum shear  stresses  are  larger  near  the  top  than  in  the  mid-depth  region, 
and  that  they  increase  again  as  depth  increases.  This  is  believed  due 
to  the  relatively  large  horizontal  stresses  which  are  induced  by  the 
deformation  mode  of  the  embankment,  (cf.  Figure  21). 

Two  contours  of  1       /7H  for  the  normal  loading  approximation  corre- 
sponding to  the  embankment  considered  are  shown  in  Figure  26  as  dashed 
lines.  They  indicate  a  shear  stress  less  than  that  produced  by  the 
elastic  embankment  in  a  shallow  zone  below  the  embankment,  but  larger 
shear  stresses  at  depth. 

Data  for  a  variety  of  embankment  shapes,  v/ith  \i   =  0.3,  and  a  =  15°, 
30°,  U50,  60°,  and  75°  are  shown  in  Figures  27  to  31,  respectively. 
In  these  figures,  the  maximum  value  of  t    &£.  &   particular  depth  is 
plotted  as  a  function  of  depth  for  various  L/H  ratios.  The  horizontal 
location  of  the  point  at  which  this  maximum  value  occurs  is  also  shown. 
It  can  be  observed  that  in  Figure  27,  as  the  L/H  ratio  increases  for 
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Figure  26.  Contours    for   Max.  Shear  Jmax/YH 
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a  =  15°,  the  peak  i         increases  in  magnitude,  and  acts  at  an  increas- 
ing depth  below  the  embankment.  The  horizontal  location  of  the  maximum 
shear  stress  at  a  particular  depth  moves  from  a  position  near  the  toe 
of  the  slope  immediately  beneath  the  embankment  to  the  center  line  of 
the  embankment  at  a  depth  which  depends  upon  the  L/H  ratio. 

A  similar  trend  is  shown  in  Figure  28  for  a  ■  30°.  However,  in 
Figure  29  (a  =  ^5°)>  the  largest  shear  stress  occurs  near  the  toe  of 
the  slope  for  L/H  =  0.0.  Although  the  smallest  L/H  ratio  shown  in 
Figures  30  and  31  is  0.5,  the  development  of  large  shear  stress  near 
the  toe  of  narrow  steep  embankments  is  clearly  indicated. 

The  dashed  line  in  Figure  29  shows  the  magnitude  of  the  maximum 
t    as  a  function  of  depth  for  the  normal  loading  approximation  corre- 
sponding to  the  i+5°  embankment  for  which  L/H  =  1.  It  is  evident  that 
the  peak  magnitudes  are  nearly  the  same  for  the  two  cases,  but  it  occurs 
at  approximately  twice  the  depth  in  the  case  of  the  normal  loading  ap- 
proximation. A  similar  effect  is  evident  in  Figure  26  for  L/H  =  3- 
Thus  the  influence  of  the  elastic  embankment  is  more  pronounced  nearer 
the  surface  where  softer  soils  might  be  expected.  As  a  result,  it  may 
be  that  current  estimates  of  stability,  potential  creep  and  other  shear 
stress  related  phenomena,  for  soils  at  shallow  depths  beneath  embank- 
ments, are  unconservative. 

Relationship  of  Results  to  In-Situ  Stresses 
It  is  not  immediately  clear  what  relationship  these  results  have 
to  stresses  which  actually  exist  in  the  field.  In  the  case  of  a  built- 
up  embankment,  it  is  likely  that  the  embankment  material  will  exhibit 
significantly  different  mechanical  properties  from  the  foundation 
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material.  For  a  cut-down  slope,  the  assumption  of  homogeneity  in  the 
two  zones  nay  be  more  nearly  justified.  The  non-linearity  in  the  me- 
chanical response  of  most  natural  materials  will  undoubtedly  also  in- 
fluence the  results.  However,  the  feature  which  may  be  most  significant, 
at  least  in  the  case  of  built-up  embankments,  is  the  fact  that  they  are 
constructed  in  layers  rather  than  instantaneously.  Thus  when  the  top- 
most lift  is  placed  on  an  earth  embankment,  the  upper  material  does  not 
undergo  strain  due  to  elastic  deformation  of  the  embankment  resulting 
from  the  stresses  imposed  by  the  entire  mass.  Rather,  the  strains  are 
due  only  to  the  increment  of  stress  imposed  by  this  layer.  The  degree 
to  which  the  results  would  be  changed  is  not  clear.  However,  it  is  be- 
lieved that  the  results  presented  herein,  provide  a  more  realistic  es- 
timate of  stress  conditions  than  that  computed  from  the  normal  loading 
approximation. 

Effect  of  Pe suits  on  Stress  Path  Determination 
Lambe  [21]  has  suggested  that  the  "stress  path"  method  for  pre- 
diction of  vertical  settlements  of  cohesive  soils  is  superior  to  con- 
ventional analyses  in  cases  where  compression  is  clearly  not  one- 
dimensional.  This  approach  involves  three  basic  steps  (Lambe,  21): 

1.  Estimation  of  the  effective  stress  path  of  an  "average" 
element  in  the  compressible  layer,  for  the  field 
loading . 

2.  Performance  of  a  laboratory  compression  test  which 
duplicates,  insofar  as  practicable,  the  field 
effective  stress  path. 

3.  Computation  of  settlement  by  multiplying  the 
thickness  of  the  layer  considered  by  the  axial 
(vertical)  strain  from  the  laboratory  test. 


Because  the  strains  in  the  laboratory  sample  depend  upon  the  applied 
stresses,  the  method  requires  a  means  of  correctly  assessing  the  in- 
situ  stresses. 

A  comparison  of  the  total  stress  paths  for  several  points  under 
the  center  line  of  an  elastic  embankment  (a  =  30°,  L/Hf.  .   =  0.5, 
u  =  0.3),  with  those  computed  using  the  normal  loading  approximation, 
is  shown  in  Figure  32.  The  dashed  "initial  stress"  line  shows  the 
state  of  stress  in  an  elastic  half-space,  for  which  u  =  0.3,  before 
construction  of  the  embankment.  The  three  points  shown  on  the  line 
correspond  to  the  stress  states  depths  of  0.5,  1.0  and  2.0  times  the 
final  height  of  the  embankment,  H„.   .. .  The  open  points  show  the 
stress  paths  during  "construction"  of  an  elastic  embankment  continuous 
with  the  foundation  material.  The  solid  points  show  the  stresses  for 
corresponding  embankment  heights,  determined  by  the  conventional  method. 
Oeveral  features  of  this  comparison  are  especially  noteworthy: 

1.  At  relatively  shallow  depths  (y/H  =  0.5),  the 
conventional  method  leads  to  a  total  stress  path 
which  lies  entirely  below  the  K  line.  That  is, 

one  would  predict  relatively  small  shear  settlements. 
However,  on  the  basis  of  the  elastic  embankment 
analysis,  the  estimated  shear  induced  settlement 
would  likely  be  larger,  and  compression  settlement 
would  be  less. 

2.  At  greater  depth  (y/H^.    =  2.0)  both  methods  lead 
to  stress  paths  which1 lie  above  the  k  line.  The 
two  paths  are  closer,  and  the  shear  stress  under 
the  elastic  embankment  is  actually  less  than  that 
due  to  the  normal  loading  approximation. 

3.  At  intermediate  depths  (y/H„.   ..  »  1.0)  the  normal 
loading  approximation  remains  relatively  close  to 
the  k  line.  The  stress  path  due  to  the  elastic 
embankment  is  still  considerably  steeper. 

Because  of  the  influence  of  the  applied  stress  path  on  the  measured 
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laboratory  settlements,  and  therefore  on  the  computed  field  settlements, 
it  would  seem  essential  to  estimate  the  predictive  capability  of  the 
method  to  the  field  stresses  accurately.  In  the  case  considered,  the 
stresses  produced  by  the  elastic  embankment  are  significantly  different 
from  those  due  to  the  normal  loading  approximation,  at  least  at  shallow 
depths.  The  effect  of  this  difference  on  the  results  predicted  by  the 
stress  path  method  is  not  obvious,  however,  this  question  would  appear 
to  deserve  further  attention. 
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PART  II 
ANALYSIS  OF  THE  VISC0ELA3TIC  EMBANKMENT 


102 


INTRODUCTION 

Many  embankments  are  composed  of  cohesive  soils.  The  deformation 
of  these  embankments  are  composed  of  three  parts:   (a)  the  immediate 
deformations  which  are  due  to  the  elastic  characteristics  of  the  em- 
bankment and  the  underlying  foundation;  (b)  long-term  deformation  which 
is  probably  attributed  to  a  creep  mechanism  in  response  to  shear  stress; 
(c)  the  deformation  due  to  consolidation  of  the  soils  underlying  the 
embankment . 

In  this  part  of  this  dissertation  a  viscoelastic  embankment  is 
analyzed.  The  material  comprising  the  embankment  and  the  foundation 
was  assumed  to  behave  in  a  linear  viscoelastic  manner. 

The  creep  curve  used  was  obtained  from  the  results  of  triaxial 
creep  tests  performed  on  kaolin  soil  from  Edgar,  Florida  compacted 
at  optimum  water  content  for  the  compaction  method  used  [30]. 

The  solution  of  the  problem  was  obtained  by  applying  the  quasi- 
elastic  approximation  proposed  by  Schapery  [32],  to  the  associated 
elastic  solution,  with  creep  function  in  place  of  elastic  constants. 
That  is,  if 

u  =  Ef (x,  y,  n) 

(122) 
v  -  Eg(x,  y,  m) 

are  the  elastic  displacements,  then  the  viscoelastic  displacements,  for 
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a  linear  viscoelastic  material,  are: 

u(T)  =  E  (T)  f[x,  y,  n(T)] 


v(T)  =  Ec(T)  g[x,  y,  n(T)] 


where : 


(123) 


E  (T)  is  the  creep  compliance  of  both  the  embankment  and  the 
foundation. 

ji(T)  is  Poisson's  ratio  in  creep. 

T  is  the  time. 

u(T)   is  the  horizontal  time  dependent  displacement. 

v(T)   is  the  vertical  time  dependent  displacement. 

If  Poisson's  ratio  in  creep  is  independent  of  time,    then  the 
quasi-elastic  method  is  exact. 
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TRIAXIAL  LOAD  AID  CREEP  TESTS 


Material 


As  mentioned  above,  the  soil  used  in  the  determination  of  the 
creep  compliance,  E  ,  is  a  moderately  well  crystallized  kaolin  soil 
from  Edgar,  Florida. 

Classification  properties  of  the  soil  are  given  in  Table  2.  The 
grain  size  distribution  and  compaction  characteristics  of  the  soil  are 
shown  in  Figures  33  and  3^ • 

Table  2 

Classification  Properties  of  the  Soil 

Liquid  Limit  58.5 

Plasticity  Index  %  22.0 

Specific  Gravity  of  Solids  2.6 

Clay  Fraction  %  78.0 

Preparation  of  Soil  Specimens 
Test  specimens  were  prepared  at  the  desired  moisture  content  by 
impact  compaction  in  a  1.3  inch  diameter  mold  with  a  special  drop 
hammer  weighing  3.62  pounds  and  height  of  fall  of  10  inches  [30J. 
Specimens  were  prepared  in  three  layers.  Impact  compaction  curves  are 
shown  in  Figure  3^.  Note  that  the  curves  are  essentially  similar  in 
appearance  to  those  of  many  naturally  occurring  silty  clays,  with  a 
moderately  well  defined  peak  point. 
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Experimental  Apparatus  and  Procedure 

The  triaxial  apparatus  used  was  of  the  "Norwegian"  type  accomo- 
dating of  1.3  inch  diameter  by  2.8  inch  long  specimen,  described  in 
detail  by  Andresen  et  al  [11.  Lateral  pressures  were  applied  through 
a  confining  fluid  to  which  pressure  is  applied  by  small  "constant 
pressure  cells".  Axial  load  is  applied  through  a  piston  by  placing 
weights  on  the  load  hanger. 

The  creep  test  was  conducted  by  permitting  the  specimen  to  come 
to  equilibrium  under  a  desired  confining  fluid  pressure.  After  the 
specimen  reached  equilibrium,  the  creep  stresses  were  applied  axially. 

The  axial  displacement  of  the  specimen  was  determined  by  measuring 
the  vertical  displacement  of  the  cross  bar  of  the  loading  frame  by  a 
linear  variable  differential  transformer  (LVDT)  which  converts  displace- 
ment into  an  electrical  signal. 

The  radial  displacement  of  the  specimen  is  determined  by  a  radial 
displacement  sensor  which  is  a  modification  of  the  lateral  strain  in- 
dicator suggested  by  Bishop  and  Henkel  [3],  and  is  shown  schematically 
in  Figure  35.  The  apparatus  and  experimental  techniques  are  discussed 
in  detail  by  Perloff  [30], 
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Figure  35.    Schematic    Diagram    of    Radial 
Displacement    Sensor    and    Spe- 
cimen  on    Triaxial    Cell    Base  (Perloff) 
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RESULTS  AND  DISCUSSION  OK  CONFINED  CREEP  TESTS 

Axial  Strain 
A  typical  creep  curve  plotted  to  logarithmic  scales  is  shown  in 
Figure  36.  This  figure  indicates  the  axial  strain,  €.  ,  as  a  function 
of  time  for  an  impact  compacted  kaolin  specimen.  This  result  is  typi- 
cal of  the  tests  which  have  been  run  [30],  including  tests  of  1000 
minutes  duration.  It  is  clear  that  the  creep  curve  appears  as  a 
straight  line  on  the  logarithmic  plot.  The  equation  of  such  a  line  is: 


log  e,  =  b  log  T  +  log  a 


(121;) 


e±   =  aT 


(125) 


where  T  is  time,  and  a  and  b  are  constants  of  the  equation.  The 
constant  a  is  readily  determined  as  the  magnitude  of  the  creep  strain 
at  a  time  of  one  minute.  The  constant  b  is  the  logarithmic  slope  of 
the  straight  line,  i.e., 


log 


log 


K> 


TCT 


(126) 


Because  the  materials  comprising  the  embankment  and  the  foundation 
are  assumed  to  behave  in  a  linear  viscoelastic  fashion,  a  and  b  must  be 
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assumed  independent  of  the  stress  level.  Perloff  [30]  has  shown  that 
the  constant  b  is  indeed  independent  of  the  stress  level,  but  the  con- 
stant a  is  not,  at  least  for  levels  exceeding  30$« 

The  ratio  of  the  axial  to  the  radial  strains  is  shown  in  Figure  37. 
This  figure  indicates  that  this  ratio,  which  is  equal  to  Poisson's  ratio 
u,  is  independent  of  tine,  i.e.,  the  relationship  between  Poisson's  ratio 
and  tine  is  a  straight  line  parallel  to  the  tine  axis.  Therefore,  the 
stresses,  which  are  functions  of  Poisson's  ratio  are  tine-independent. 


Creep  Compliance 
The  creep  compliance  can  be  determined  from: 


a     -  a_ 
E  =  -± £  (127) 


Combining  equations  (125)  and  (127) 


o  -  a 
E  =  -i 3  (123) 

aT" 


To  determine  the  constant  a  and  b,  consider  the  correspondence 
between  e  and  T  in  equation  (125). 

-k 

1  -  at  T  =  1  min.,    e  =  9.2h  x  10 

Cubstituting  into  equation  (125): 
a  =  9.2^  x  10 

2  -  at  T  «=  10  min.,   e  =  9.8l  x  10_  +  .  From  equation  (125): 

b  =  0.0259 
Therefore : 
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a,  -  a, 


Ec  =  " 1     0  025Q  ^129) 


Let  E   be  the  creep  compliance  at   T  =  0.001  rain.  Substituting  into 
o 


equation  (129) 


gl  -  g3 


s  =  -tf ,  (130) 


c 


ZX 


~o    9.2U  x  10   x  0.8358 
Dividing  equation  (12Q)  over  equation  (130),  one  obtains: 


Ec        _       0.8353  ,__., 

~r~  -    To.o259  ^  °  ; 

o 


where  T  has  to  be  in  rainutes, 


is-  =  ^8 ,   } 

^T      11112  t0-0259  < 32; 

o 


where  T  has  to  be  in  hours. 
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RESULTS  AMD  DISCUSSION 


Introduction 


The  displacement  within  and  under  a  long  viscoelastic,  isotropic 
and  homogeneous  embankment  continuous  with  the  underlying  material  of 
similar  properties,  has  been  investigated  for  a  =  30°,  l/H  =  1,  u  =  0.3 
and  creep  characteristics  shown  in  Figure  36. 

The  magnitude  of  Poisson's  ratio  was  chosen  to  be  0.3,  which  is 
different  from  the  test  data  given  in  Figure  37.  This  discrepancy 
occurs  because  the  test  was  conducted  on  a  compacted  unsaturated  speci- 
men, while  the  material  underneath  the  embankment  is  likely  to  be  satu- 
rated. The  value  u  =  0.3  was  considered  to  be  an  appropriate  average. 
Further,  it  is  assumed  that  the  embankment  is  constructed  instantaneous- 
ly. 

Effect  of  Time  on  the  Embankment  Shape 
A  typical  result  of  the  analysis  is  illustrated  in  Figure  33.  This 
figure  shows  the  change  of  the  shape  of  the  embankment  with  respect  to 

time,    in  dimensionless  form,  VE/ylf"   (where  V  =  u  +  iv).     The  dashed 

o 
lines  show  the  displacements  of  the  base  and  under  the  edge  of  the  em- 
bankment. These  results  indicate  that,  for  the  soil  considered,  the 
immediate  displacement  is  almost  sixty  per  cent  of  the  displacement 
occurring  within  the  first  ten  years.  Ninety  per  cent  of  this  displace- 
ment occurs  within  the  first  month.  The  maximum  value  of  the  displacement 


115 


,H^/°°3   x   »uaui33D|ds!a 


116 


occurs  along  the  center  line  and  decreases  as  x  H  increases  and  becomes 
very  small  at  large  x/H  values,  i.e.,  less  than  0.012  for  x/H  =  50. 

Analysis  of  Vertical  and  Horizontal  Displacement 

Figure  39  shows  the  effect  of  time  on  the  vertical  displacements 
of  five  points  from  the  embankment,  shown  schematically  in  the  figure. 
These  results  are  typical  of  all  the  points  in  the  region.  The  observed 
straight  line  logarithmic  plot  of  displacement  as  a  function  of  time  is 
implied  by  equations  (123)  and  the  experimental  data,  Figure  36. 

Figure  HO  shows  the  effect  of  time  on  the  compression  of  the  em- 
bankment discussed  above,  between  points  A  and  E.  This  figure  indicates 
that  the  vertical  displacements  induced  within  the  embankment  at  any 
time  are  approximately  fifteen  per  cent  of  the  total  displacement. 

Figure  kl   shows  the  relative  vertical  displacement  between  points 
E  and  C,  that  is,  Lie  difference  between  curves  E  and  C  in  Figure  39« 
Thus  the  base  of  the  embankment  is  settling  at  a  faster  rate  in  the 
center  than  at  the  toe. 

The  effect  of  tine  or;  the  relative  vertical  displacement  between 

the  center  top  and  the  top  edge  of  the  embankment  is  shown  in  Figure  k2. 

This  indicates  that  very  little  differential  settlement  between  points 

A  and  B,  i.e.,  for  L  =  II  =  30  ft.,  a  =  30°,  y  -   120  pounds  per  cubic 

foot  and  E  =  E  ,  the  computed  differential  settlement  between  A  and 
c    c  * 
o 

B  during  the  first  ten  years  is  0.059  inches. 

The  relative  horizontal  displacements  are  also  small,  as  indicated 
in  Figure  U3 •  This  figure  shows  the  effect  of  time  on  the  horizontal 
displacements  of  five  points  from  the  embankment.  It  is  clear  that  the 
horizontal  displacement  at  points  A  and  E  is  zero,  this  is  due  to  the 
symmetry  of  the  problem. 
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PART  III 
APPLICATION  OF  METHOD  TO  ANALYSIS  OF  AN  EXCAVATION 
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IMTRODUCTION 

The  method  presented  above  can  solve  many  plane  strain  problems  of 
homogeneous,  isotropic  material  and  can  be  extended  to  include  aniso- 
tropic material  and  layered  systems.  In  fact,  the  method  does  not  dif- 
ferentiate in  any  way  between  embankments  and  foundations,  or  between 
slopes  and  cuts.  However,  the  transformation  formula,  which  transforms 
the  geometric  boundary  into  a  half-space  differs  in  each  case. 

As  an  illustration  of  the  application  of  the  method,  the  stresses 
induced  within  an  elastic,  homogeneous  and  isotropic  mass  due  to  exca- 
vation of  an  open  trench  in  a  semi-infinite  mass,  and  the  application 
of  a  uniform  strip  load  on  the  base  of  the  excavation  are  analyzed. 

The  problem,  illustrated  in  Figure  kk   (a),  is  that  of  a  long  com- 
pensated foundation,  in  which  a  depth  of  soil  is  excavated  to  reduce  the 
net  load  applied  to  compressible  subsurface  materials.  The  usual  method 
for  estimating  the  stress  distribution  is  to  assume  that  the  net  load 
applied  at  the  base  of  the  excavation  is  equal  to  the  foundation  load, 
less  the  unit  weight  of  soil  times  the  depth  of  excavation,  and  that 
this  stress  distributes  in  the  same  fashion  as  a  strip  load  on  the  sur- 
face of  a  half-space  placed  at  the  elevation  of  the  base  of  the  excava- 
tion. This  is  referred  to  herein  as  the  Boussinesq  approximation,  and 
is  illustrated  in  Figure  hk   (b). 
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THE  TRANSFORMATION  OF  THE  GEOMETRIC 
BOUNDARY  INTO  A  HALF- SPACE 


Applying  the  Schwarz-Christoffel  transformation  (56)  to  transform 
Figure  k3   (a)  into  Figure  k5   (b),  leads  to: 


z     L_ 
H  =  HP 
0 


t  _  1 


/   A  -  p   .N  ,,  .  . 


(133) 


where: 


p  is  the  modulus  of  the  integral. 
F  is  the  value  of  the  integral  at  point  p. 
Figures  k6   and  1+7  give  the  value  of  p  as  a  function  of  L/H  or  H/L. 
Knowing  p,  the  value  of  F  is  determined  from  Figure  U8. 
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RESULTS  AND  DISCUSSION 

Figure  1+9  shows  the  reduction  in  vertical  stress,  a   /7H,  along  the 
center  line  due  to  the  excavation,  calculated  by  this  method  and  by  the 
Eoussinesq  approximation  for  L/H  ■  1  and  u  =  O.U.  Here  again,  the  stress 
determined  by  the  Boussinesq  analysis  is  independent  of  Poisson's  ratio; 
the  stress  determined  herein  depends  upon  u.  Figure  U9  indicates  that 
the  reduction  in  vertical  stress,  at  depths  between  y/H  =  0  and  y/H  =  1, 
is  less  than  that  computed  by  the  Boussinesq  approximation  and  higher 
elsewhere.  The  difference  between  the  computed  vertical  stress  reduc- 
tion and  Boussinesq  value  along  the  center  line,  expressed  as  a  percen- 
tage of  the  Boussinesq  value,  is  shown  in  Figure  50,  for  L/H  =  0.5,  1 
and  2.  For  all  L/H  ratios,  the  reduction  in  vertical  stress,  between 
the  bottom  of  the  foundation  (y/H  =  0)  and  approximately  y/H  =  1,  is 
slightly  smaller  than  computed  by  the  normal  loading  approximation,  but 
is  significantly  larger  below  the  depth  y/H  =1.  Furthermore,  this 
figure  shows  that  l/h  ratio  has  pronounced  effect  on  the  difference 
between  both  methods;  a  smaller  L/H  ratio  leads  to  a  greater  disparity 
between  the  results  of  the  two  methods. 

The  reduction  in  vertical  stress  under  the  edge  of  the  excavation, 
is  shown  in  Figure  51.  The  vertical  stress  reduction  is  significantly 
higher  than  that  computed  by  the  Boussinesq  approximation.  The  dif- 
ference becomes  especially  apparent  at  the  bottom  of  the  foundation, 
where  the  vertical  stress  reduction  is  approximately  1.8  times  that 


131 


H      I 


H 


%$. 


308Z&0 


I      x/H 
fy/H 

Excavation  Analysis 


r 


n 


T&z&j&gsam 


WA-W- 


q-yH 


,  i 


x/H 


|y/H 

Boussinesq   Approximation 


x 


2- 


°    * 

o 
u 

X 

UJ 


4- 


irt     5 
O 
ffi 

o   6 
m 

f    7 

a. 
a> 
O 

8 


Reduction  in  Stress,  Ap(ay/yH) 
0.1        0.2      0.3      0.4      0.5       0.6      0.7       0.8      0.9      1.0 


1          1          1 
^■0.4 

1             1             1            1             i 

1         > 

- 

L/H-I.O 

- 

— 

/    /     ^^Excavation 

- 

- 

\Y                   Analysis 
'     \v_  Boussinesq 

- 

- 

f 

Approximation 

- 

- 

i    /    /  i            i 

i           I           i           i            i 

i 

10 

Figure  49.   Reduction  in  Vertical    Stress   Along  Center 
Line    Due    to    Excavation 


132 


-10 
0 


Net    Reduction  (%   of    Boussinesq    Value) 

0  10  20         30        40         50         60  70         80 


Figure  50.     Net    Reduction    in    Computed  Vertical   Stress 
Along   Center   Line    Due  to    Excavation, 
Compared   to   Boussinesq   Value 


133 


uLlL   .i 


H       I 


^w^1- 


x/H 

fy/H 

Excavation   Analysis 


L 

L 

"1" 

1 

IT 

i 

i 

I 

l 

q  =  /H 

VaV^\VvWw>v<W 

WWi 

^W 

x/H 
♦  y/H 

Boussinesq    Approximation 


Reduction    in    Stress,    Ap(<ry//H) 
O.I         0.2       0.3        0.4       0.5       0.6       0.7       0.8       0.9 


.0 


Figure  51.   Reduction    in    Vertical    Stress    Under    Edge 
of    Excavation 


134 


usuaHy  assumed. 

Figure  52  shows  the  percentage  differences  between  the  excavation 
analysis  and  Boussinesq  value,  below  the  edge  of  the  excavation,  for 
L/H  a   0.5,  1  and  2,  and  a  =  90°.  It  can  be  seen  clearly  from  this  figure 
that  the  difference  is  larger  at  the  bottom  of  the  foundation  and  de- 
creases as  y/H  increases  reaching  a  minimum  value  at  approximately 
y/H  =  1.  At  greater  depth,  this  figure  shows  only  a  slight  increase 
in  the  difference  for  high  value  of  L/H  and  significant  increases  in 
the  difference  for  small  value  of  L/H,  and  the  difference  is  almost 
constant  below  y/H  =  10. 

Figure  53  shows  contours  of  vertical  normal  stress,  o\  /q,  due  to  a 
uniform  strip  load,  q,  placed  on  the  base  of  the  foundation,  for  L/H  =  1, 

u  =  O.k   and  q  =  1.  The  dashed  lines  show  stress  contours  in  terms  of 
air/q,  for  t..e   Boussinesq  approximation.  It  is  clear  that  the  stress  is 
smaller  in  the  vicinity  of  the  loaded  area  than  that  computed  by 
Boussinesq  value  an J  larger  elsewhere,  this  is  a  reverse  picture  of  the 
stress  reduction  shown  in  Figures  k$   to  52. 

The  net  increase  in  vertical  normal  stress  due  to  the  combined 
effect  of  excavation  and  application  of  stress  at  the  base  of  the 
compensated  foundation  is  shown  in  Figure  5^>  for  L/H  =  1,  u  =  O.k 
and  q  =  l.UyH.  The  results  depend,  of  course,  upon  the  magnitude  of  q 
relative  to  the  depth  of  the  excavation.  A  value  of  q  =  l.kyll   is 
representative  of  that  for  a  multistory  structure.  Figure  5k   (a)  shows 
that  the  net  increase  in  vertical  stress  along  the  center  line  has  the 
same  value  at  y/H  =  0,  but  it  dissipates  faster  than  that  computed  by 
the  Boussinesq  analysis.  Figure  5k   (b)  shows  that  the  net  increase  in 
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vertical  stress  under  the  edge  of  the  excavation  is  smaller  than  that 
usually  assumed  especially  at  the  bottom  of  the  foundation  where  the 
Boussinesq  value  is  35  per  cent  greater  than  the  calculated  value. 
Thus,  one-dimensional  settlements  calculated  using  the  Boussinesq  ap- 
proximation will  be  larger  than  those  determined  by  the  excavation 
analysis. 
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CONCLUSION  AND  PEC  OMKEI  DAT  IONS 


Conclusions 


The  conclusions  to  be  drawn  from  the  analysis  presented  herein  fall 
into  three  major  categories:  those  relating  to  the  method,  those  re- 
lating to  the  elastic  results  and  those  relating  to  the  viscoelastic 
results. 

The  conclusions  drawn  from  the  method  are : 

1  -  The  method  developed  permits  solution  of  a  wide  variety  of  plane 

strain  problems  including: 

a)  Plane  strain  problems  with  free  boundary  loading. 

b)  Plane  strain  problems  with  external  static  stresses  applied 
to  the  boundary. 

c)  Plane  strain  problems  where  the  displacements  of  points  of  the 
boundary  is  given. 

d)  Mixed  plane  strain  problems. 

e)  Plane  strain  problems .  where  its  boundary  is  under  dynamic  effect, 

2  -  In  order  to  apply  this  method,  the  stresses  have  to  be  defined  as 

a  result  of: 

a)  The  stress  due  to  the  geometric  shape  of  the  boundary, 
b).  The  stress  due  to  the  weight  of  the  medium. 

The  conclusions  drawn  from  the  results  of  the  elastic  analysis  are: 
1  -  The  distribution  of  stresses  within  and  beneath  an  elastic  embank- 
ment. 


lto 


a)  The  horizontal  distribution  of  vertical  stress  beneath  the 
embankment  is  more  nearly  uniform  than  is  usually  assumed. 
Thus,  the  differential  settlements  computed  using  the  normal 
loading  approximation  will  be  larger  than  those  determined  using 
the  stress  distribution  presented  herein. 

b)  The  horizontal  and  vertical  shear  stresses  created  in  the  foun- 
dation material  by  the  embankment,  are  found  to  significantly 
higher  at  shallow  depths  for  the  elastic  embankment  than  for 
the  normal  loading  approximation. 

c)  The  horizontal  location  of  the  maximum  shear  stress  at  a  par- 
ticular depth  moves  from  a  position  near  the  toe  of  the  slope 
immediately  beneath  the  embankment  to  the  center  line  of  the 
embankment  at  a  lepth  which  depends  upon  the  L/H  ratio. 

2  -  The  distribution  of  stresses  within  an  elastic  foundation. 

a)  The  reduction  in  the  vertical  normal  stress  due  to  the  excava- 
tion is,  in  general,  larger  than  that  computed  by  the 
Boussinesq  analysis. 

b)  The  net  increase  in  vertical  normal  stress  due  to  the  combined 
effects  of  excavation  and  application  of  stress  at  the  base  of 
a  compensated  foundation  is  smaller  than  that  usually  assumed. 
Thus,  the  one-dimensional  settlement  calculated  by  this  method 
will  be  smaller  than  those  determined  by  the  Boussinesq  analy- 
sis. 

The  conclusions  drawn  from  the  viscoelastic  results,  for  the  soil 
tested,  indicate  that  sixty  per  cent  of  the  total  displacement  which 
occurs  in  the  first  ten  years  is  due  to  the  elastic  effect,  however, 
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ninety  per  cent  of  this  displacement  occurs  in  the  first  month,  and  most 
of  it  is  in  the  foundation. 

Recommendations 

1  -  The  method  presented  herein  can  be  extended  to  include : 

a)  The  effect  of  a  difference  in  the  stiffness  of  the  embankment 
and  the  foundation  on  the  distribution  of  stresses  and  dis- 
placements. 

b)  An  anisotropic  embankment  resting  on  an  anisotropic  foundation 
material. 

c)  The  layered  system  embankment  resting  on  layered  foundation 
system. 

d)  Tr.e   effect  of  a  dynamic  load  on  the  boundary  of  the  embankment, 

e)  The  non-linear  viscoelastic  embankment. 

2  -  The  application  of  the  method  presented  herein  on  different  plane 

strain  problems  is  highly  recomraende  ; . 

3  -  The  possibility  of  verifying  the  results  included  in  this  investi- 

gation, experimentally,  should  be  studied. 
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APPENDIX  A 

Typical  Test  Data 
Soil  -  kaolin 
Initial  height  (L  )  =  2  i||  in 


Diameter  (D  )       =1  ^pF1  i° 


Compaction  -  impact  -  3  layers  -  15  blows  each  layer 
Confining  pressure  (a. )   =  3  kg  cm 

Existing  axial  load     =  0 
(Area)(a1  -  a3)Q 

Additional  axial  load   =18.1+ 
(Area)  (a,  -  aj 

Initial  net  weight  of  the  specimen  =  119.32 
Final  net  weight  of  the  specimen  =  H9J+6 
Dry  weight  of  the  specimen  «=  Qk.06 
Initial  moisture  content  =  27t,j 
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APPENDIX  B 


COMPUTER  PROGRAM  FOR  CALCULATION  OF 
THE  TRANSFORMATION  INTEGRAL  CONSTANTS 


$ID   338U*2*10**BALADI*GEORGE 

$EXECUTE       IBJOB 

$IBJOB        GO 

ilBFTC  TRIAL 

C     M  MUST  BE  ODD 

COMPLEX  LF,F,F1,FM,EI(20) ,D1,D2,EIF,EIG,RIA,EIB,FF,CIBC 

REAL  N 
200  READ(5,99)  DBET,NN 
99  FORMAT (Flo. 5 ,110) 

WRITE (6, 11) 

11  F0RMAT(15X,3HEIF,25X,3HEIG,25X,3HEIA,22X,3HEIB) 

100  READ(5,101)  N 

101  FORMAT (F10. 5) 
WRITE (6, 12)  N 

12  FORMAT (5X,2HN=,F10. 5) 
READ(p,Ul)  BETT 

1*1  F0RMAT(F10.5) 
DO  109  J=1,NN 
PP=J-1 

BET=BETT+PP*DBET 
WRITE (6, 13)  BET 

13  FORMAT (5X,5HBETA=,F10. 5) 
K=l 

M=251 
T1=0.0 
DELT=0 . 9*BET 

102  D1=(0. 0,0.0) 
D2=(0. 0,0.0) 
DO  104  1=1, M 
P=M-1 

QfI-1 

T=T1+DELT*Q/P 

FF= ( T**2 - 1 . 0 ) / ( T**2 -BET**2 ) 

if(real(ff) .eq.o.o.and.aimag(ff) .eq.o.0)f=(o. 0,0.0) 
if(real(ff).eq.o.o.and.almag(ff).eq.o.o)go  to  190 
lf=clog(ff)/n 

F=CEXP(LF) 
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190  IF(I.EQ.I)   F1=F 

IF(I.EQ.I)   GO  TO  104 
IF(I.EQ.M)  FMcF 
IF(l.EQ.M)   GO  TO  104 
IF((l-l)/2.EQ.l/2)   GO  TO  103 
D2=D2+F 
GO  TO  104 

103  M=D1+F 

104  CONTINUE 

EI ( K ) > ( DELT/ ( 3 . 0*P ) ) * ( Fl+2 .  C^Dl+U . 0*D2+FM ) 
K=K+1 

GO  TO  (105,105,106,111,121,122,126,123,124,112,125,113,114,116,110 
1),K 

105  T1=0.9*BET 
DELT=0.09*BET 
GO  TO  102 

106  T1=0.99*BET 
DELT=0.009*BET 
GO  TO  102 

111  T1=0.999*BET 
DELT=6.0009*BET 
GO  TO  102 

121  T1=0.9999*BET 
DELT=0.00009*BET 
GO  TO  102 

122  T1=0.99999*BET 
DELT=0 . 000009*BET 
GO  TO  102 

126  T1=0.999999*BET 

DELT=0 . 0000009*BET 
GO  TO  102 

123  T1=1.0000001*BET 
DELT=0 . 0000009*BET 
GO  TO  102 

124  T1=1.000001*BET 
DELT=0 . 000009*BET 
GO  TO  102 

112  T1=1.00001*BET 
DELT=0.00009*BET 
GO  TO  102 

125  T1=1.0001*BET 
DELT=0.0009*BET 
GO  TO  102 

113  T1=1.001*BET 
DELT=0.009*3ET 
GO  TO  102 

114  T1=1.01*BET 
DELT=0.09*BET 
GO  TO  102 

116  T1=1.10*BET 

DELT=1.0-1.1*BET 
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GO  TO  102 

110  eif=ei(1)+ei(2)+ei(3)+ei(u)+ei(5)+ei(6)+ei(7) 

eig=eif+ei(8)+ei(9)+ei(10)+ei(11)+ei(12)+ei(13)+ei(1u) 

eia=eif/eig 

eib=eig/eif 

cib=aimag(eig) 

cibc=eif/cib 
109  wkite(6,108)  eif,eig,eia,eib,cibc 
108  format(5x,2f10.5,5x,2f10.5,5x,2f10.5,5x,2f10.5,5x,2f10.5) 

GO  TO  200 

END 
$DATA 
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APPENDIX  C 


COMPUTER  PROGRAM  FOR  THE  DETERMINATION  OF 
STRESSES  AND  DISPLACEMENTS  WITHIN  AND  UNDER 
LINEAR  ELASTIC  AND  VISCOELASTIC  EMBANKMENTS 


$LD  338U*9*50**BALADI*GEORGE 

$EXECUTE  IBJOB 

$IBJOB  GO 
$IBFTC  BALADI 

C  LIST  OF  SYMBOLS 

C  KXKXXXKXOOCXXXXXXKMXXKKKMXXXKXXXKXMXKMXXKXHM 

C  *  PO=SQRT(-1.0)                         * 

C     * * 

C  *  ZZ(II)=MU  IS  THE  POISSON'S  RATIO         * 

C     *-- * 

C  *  CC=H/L,   WHERE  H  IS  THE  HEIGHT  OF  THE            * 

C  *  EMBANKMENT  AND  2L  IS  THE  WIDTH  OF  THE         * 

C  *  TOP  OF  THE  EMBANKMENT                                               * 

C  * * 

C  *  Nl  IS  THE  ANGLE  BETWEEN  THE   SIDE  OF  THE     * 

C  *  EMBANKMENT  AND  THE  X-AXIS                                      * 

C  * * 

C  *  A(J,K)  =  (VERTICAL  STRESS)/ (GAMA*ii)                    * 

C  * * 

C  *  B(J,K)=(3iEAR   STRESS)/ (GAMA*H)                           * 

C  * * 

C  *  S( J, K)=( HORIZONTAL  STRESS) /(GAMA*H)               * 

C  * * 

C  *  BMAX(J,K)=(MAX.SHEAR)/(GAMA*H)                           * 

c  * * 

C  *  AMAX(J,K)«(MAX.oTRESS)/(GAMA*H)                        * 

c  * * 

C  *  SMIN(J,K)=(MIN.STRESS)/(GAMA*H)                        * 

c  * * 

C  *  ALFA1   IS  THE  ORIENTATION  ANGLE  OF  MAJOR     * 

C  *  PRINCIPAL  STRESS  MEASURED  COUNTER-                * 

C  *  CLOCKWISE  FROM  THE  HORIZONTAL                           * 

C  * * 

C  *  ALFA2  IS  THE  ORIENTATION  ANGLE  OF  MINOR  * 

C  *  PRINCIPAL  STRESS  MEASURED  COUNTER-       * 

C  *  CLOCKWISE  FROM  THE  HORIZONTAL           * 

C     * * 
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C  *  V(J,K)=(VERT.DISPIACEMENT*E)/(GAMA*H**2)   * 

C  * * 

C  *  U(J,K)=(H0RIZ.DISPLACEMENT*E)/(GAMA*H**2)* 

C  * * 

C  *  E  IS  THE  ELASTIC  MODULUS  * 

C  * * 

C  *  EC  IS  THE  CREEP  COMPLIANCE  * 

Q  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXEXXXXX 

COMPLEX  P0,LAM(15,20),EI(15,20),EIC(15,20),DLJ,DLK,PHIE(15,20),DEL 
1T(15,20),SPHI(15,20),RI(15,20),SHI(15,20),GI(15,20) 
REAL  N,N1,ZZ(U) 

COMMON/TRA/NN,MM,DLJ,DLK,N,C,BETA,CC,LAM,EI,EIC,PO,PR,DEJ,DEK 
COMMON/FUN/PHIE ,DELT 

COMMON/STR/A(l5,20)fB(l5,20),S(l5,20),BMAX(l5,20),AMAX(l5,20),SMIN 
1(15,20) ,ALFA1(15, 20), ALFA2(15, 20) 
COMMON/PH/SPHI , SHI 
COMMON/GES/GI 

C0MM0N/DIS/U(15,20),V(15,20) 
COMMON/POIS/ZZ 
COMMON/VET/NL ,DT , DET 
COMMON/CHIF/ADJ 

READ(5,11)  LL 

11  FORMAT ( 110) 

READ(5,12)  (ZZ(II),II«1,LL) 

12  FORMAT (2F10. 5) 
READ(5,1^)  NL,DT,DET 

Ik   F0RI«JAT(I10.2F10.5) 

100  READ(5,101)  NW,MM,DEJ,DEK 

101  FORMAT ( 2110, 2F10. 5) 
READ(5,102)  N  C,BETA,CC 

102  FORMAT (UF10. 5) 
READ(5,103)  PR 

103  FORMAT (F10. 5) 
READ(5,l6)  ADJ 

16  FORMAT (F10. 5) 
Nl=l80.0/N 
P0=1.0*(0. 0,1.0) 

CALL  TRANSF 

C     STRESS  CALCULATION 
C     

DO  13  11=1, LL 

C     THE  STRESS  FUNCTION 
C     - 


CALL  FUNCT(II) 
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CALL  STRESS (II) 

C     THE  INTEGRAL  OF  THE  STRESS  FUNCTIONS  FOR  THE  DISPLACEMENTS 
C     

CALL  PM(II) 

CALL  PSI(II) 

C     THE  ELASTIC  DISPLACEMENT 
C     

CALL  DISPLA(II) 

C     THE  VTSCOELASTIC  DISPLACEMENT 
C     

CALL  VISCO(II) 

WRITE(6,21)  NN,MM,DEJ,DEK 

21  F0RMAT(///1H0,3HNN=,I3,9X,3HMKI=,I3,9X,UHDEJ=,F9.5,UX,UHDEK=,F9.5) 
WRITE (6, 22) 

22  FORMAT (LH1,5X,  10*+ HSTRESS  DISTRIBUTION  WITHIN  AND  UNDER  LONG  ELASTI 
1C  AND  VISCOELASTIC  EMBANKMENT  WITH  FREE  BOUNDARY  LOADING) 

WRITE(6,23)  ZZ(II) 

23  FORMAT (///10X,3HMU=,F10. 5) 
WRITE(6,2U)  N1,C,BETA,CC 

2U  F0RMAT(5X,3HN1=,F10.5,5X,2HC=,F10.5,5X,5HBETA=,F10.5,5X,UHH/L=,F10 

1.5) 

WRITE (6, 25)   ((EI(J,K),A(J,K),B(J,K),S(J,K),AMAX(J  K) ,BMAX( J,K) ,SM 
1IN(J,K),ALFA1(J,K),ALFA2(J,K),U(J,k),V(J,K),K=1,MM),J=2,NN) 
25  FGira-IAT(6X,3HX/H,7X,3Hy/H,5X,8}iSIGMA(Y),3X,5IISIJEAR,3X,8HSIGMA(X),lX 
1,10!ISIGMA(M/«),1X,9HMAX.SHE^,LX,10^ 
2A2,5X,UHU(X),6X,1+HV(Y)///(LX,10F10.5,2F9.U)) 
13  CONTINUE 
GO  TO  100 
END 


$IBFTC   SUB1 

SUBROUTINE  TRANSF 

C            THE  TRANSFORMATION  FROM  THE  Z-PLANE  ONTO  THE  T-PLANE 
C  - 

COMPLEX  P0,LAM(15,20),EI(15,20),EIC(15,20),DLJ,DLK,D1,D2,ED1,ED2,T 
1,FF,LF,F,EF,F1,EF1,FM,EFM 
REAL  N,N1 

C        XXHXKKXXXXXXXXXXXKKKMKXMXXXMXKXXKXXXMMXXXXMXMMKXXXXXXKXXXXXKMXXKKX 

C     *  BETA  IS  THE  MODULUS  OF  THE  EMBANKMENT  TRANSFORMATION  INTEGRAL  * 
C     * * 
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C  *  EI(J,K)=X/H+P0*Y/H  THE  TRANSFORMATION  FORMULAE               * 

C     * * 

C  *  EIC(J.K)  IS  THE  HARMONIC  CONJUGATE  OF  El(j,K)  AT  THE  BOUNDARY  * 

C     * * 

C  *  1/C,  IS  A  CONSTANT  PRECEDE  THE  TRANSFORMATION  INTEGRAL  El(j,K)  * 

Q  XXXKKHXXXXKMKKXXMXMXXXKXXXXXKXXKKMKyMKKXKMXXKKXMKKXMKXXXXXKXKKMXXX 

C            M  MUST  BE  ODD 
C  

C0MMON/TRA/NN,l-M,DLJ,DLK,N,C,BETA,CC,LAM,EI,EIC,PO,PR,DEJ,DEK 

LAM(1,1)=(0. 0,0.0) 

EI(1,1)=-P0 

EIC(1,1)=P0 

DLJ=DEJ*PO 

DLK=-DEK*(P0**2) 

DO  10U  J=2,NN 

PP-J-1 

DO  10U  K=1,MM 

W=K-1 

LAM  ( J ,  K )  =PR*DU+  ( PP- 1 . 0 )  *DLJ+W*DLK 

IF(J.EQ.2.AND.K.EQ.l)   GO  TO  105 

M=101 

GO  TO  106 

105  G1=100.0*PR 
M=2*IFLX(G1)+1 

106  P=M-1 
D1«(0. 0,0.0) 
D2=(0. 0,0.0) 
ED1=(0. 0,0.0) 
ED2=(0. 0,0.0) 
DO  107  1=1, M 
Q=i-1 

IF(K.NE.l)  GO  TO  108 
IF(J.EQ.2)  GO  TO  109 
T=LAM(J-1,1)+DLJ*^/P 
GO  TO  110 

109  T=LAM(l,l)+PR*DLJ*Q/P 
GO  TO  110 

108  T=LAM(J,K-l)+DLK*Q/P 

110  FF=(T**2-1.0)/(T**2-BETA**2) 
LF=CLOG(FF)/N 

F=CEXP(LF) 

EF=REAL ( F ) -AIMAG ( F ) *P0 

IF(I.EQ.I)  F1=F 

IF(I.EQ.I)   EF1=EF 

iWl.EQ.l)   GO  TO  107 

IF(I.EQ.M)   FM=F 

IF(I.EQ.M)   EFM=EF 

LF(I.EQ.M)   GO  TO  107 

IF((l-l)/2    .EQ.    1/2)   GO  TO  111 
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D2=D2+F 
ED2=ED2+EF 
GO  TO  107 

111  D1=D1+F 
ED1=ED1+EF 

107  CONTINUE 

IF(K.EQ.l)  GO  TO  112 

ei(j,k)=(dlk/(3.0*p))*(f1+2.0*d1+u.0*d2+fm) 

ei(j,k)=(1.0/(cc*c))*ei(j,k)+ei(j,k-1)+p0 

ei(j,k)=ei(j,k)-po 

EIC(J,K)=(DLK/(3.0*P))*(EFl+2.0*EDl+i4.0*ED2+EFM) 

eic(j,k)=(1.0/(cc*c))*eic(j,k)+eic(j,k-1)-p0 
eic(j,k)=eic(j,k)+po 

GO  TO  10U 

112  if(j.eq.2)  go  to  113 
ei(j,1)=(dlj/(3.0*p))*(f1+2.0*d1+14.0*d2+fm) 
ei(j,i)=(i.o/(cc*c))*e:i(j,i)+ei(j-i,i)+po 

EI(J,1)=EI(J,1)-P0 

EIC(J,l)=(DLJ/(3-0*P))*(EFl+2.0*EDl-^.0*ED2+EFM) 

EIC(J,lWl.O/(CC*C))*EIC(j,l)+EIC(«J-l,l)-PO 

EIC(j,l)=£IC(J,l)+PO 

GO  TO  10U 

113  El(2,K)=(PR*DLJ/(3.O*P))*(Fl+2.0*Dl+i4.0*D2+FM) 
EI(2,K)=(1.0/(CC*C))*EI(2,K)+EI(1,1)+P0 
EI(2,K)=EI(2,K)-P0 

EIC  ( 2 ,  K )  =  ( PR*DLJ/  ( 3 .  0*P  )  )  *  ( EF1+2 .  0*EDHJ| .  0*ED2+EFM  ) 
EIC(2,K)=(1.0/(CC*C))*EIC(2,K)+EIC(1,1)-P0 
EIC(2,K)=EIC(2,K)+P0 
10*4  CONTINUE 
RETURN 
EI© 

$IBFTC  SUB2 

SUBROUTINE  FUNCT(ll) 

C            THE  CALCULATION  OF  THE  STRESS  FUNCTIONS 
C  

COMPLEX  P0,LAM(15,20),EI(15,20),EIC(15,20),DLJ,DLK,D1,D2,ED1,ED2,T 
1,FF,LF,F,EF,F1,EF1,FM,SFM,CD,DCD,CE,CCE,PEN,BCE,DD,EE,R,RR,PHI1,PH 
2I2,PHI3,PHI9)PHIi4,Pm5,PHl6,PHI7,Pra8,PHIE(l5,20),DELTl,DELT2,DELT 
33,r»ELTi4,DELT5,DELT6,DELT7,DELT8,DELT9,DELTll,DELT12,DELT(l5,20) 

REAL  N,Nl,ZZ(lt) 

COMM0N/TRA/NN,MM,DLJ,DLK,N,C,BETA,CC,LAM,EI,EIC,P0,PR,DEJ,DEK 

COMMON/ FUN/ PHIE , DELT 

COMMON/POIS/ZZ 

F0U=ZZ(II)/(1.0-ZZ(II)) 

DO  201  J=2,NN 

PP=J-1 
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DO  201  K=1,MM 

W=K-1 

CD-(LAM(J,K)««2-1.0)/(LAM(J,K)**2-BETA«*2) 

DCD=CLOG(CD)/N 

CE=(l.O/(CC*C))*(CEXP(DCT>)) 

cce=real ( ce ) -ab4ag ( ce ) *p0 
pen=real(ei(j,k))-aimag(ei(j,k))*po 

BC=1.0/((2.0*3.1'+1593)*(l.O-BETA)) 
BCB=BC*PO 

dd=cl0g((lam(j,k)+beta)/(lam(j,k)+1.0)) 
ee=clog((lam(j,k)-i.o)/(lam(j,k)-beta)) 
r=(i.o-beta)/((lam(j,k)+i.o)*(lam(j,k)+beta)) 
rr=(1.0-beta)/((lam(j,k)-1.0)*(lam(j,k)-beta)) 

c    phi  function 

C     - 

IF(N.LE.U.O)  GO  TO  202 

PHIl=2.0*3.1'+1593*(l.O/N)*PO 

PHI2»2.0*3.lUl593*(-1.0/N)*PO 

PHI3=CEXP(PHI1) 

PKE9=CEXP(PHI2) 

GO  TO  203 

202  PHI1=2. 0*3. 141593 *(-1.0/N)*P0 
PHI2=2 .0*3 .141593* (l .0/N)*PO 
PHI3— CEXP(PHIl) 
PHI9=-CEXP(PHI2) 

203  PHIU=(-BCB/2.0)*(DD+LAM(J,K)*DD-BETA+1.0) 
PHI5=(-BCB/2.0)*(EE-LAM(J,K)*EE+BETA-1.0) 

PHl6=:(-P0/(U.0*3.1i+1593))*<:L0G((LAM(j,K)-BETA)/(LAM(j,K)+BETA)) 
PHI7= ( PHlU*PHI3+PHI5*PHI9+PHl6 ) * ( 1 . O-FOU ) 

PHI8=  ( PHlit+PHI5+PHl6 )  *  ( 1 . 0+FOU  ) 

PKLE ( J , K ) =PHI7+PHI8 

C            PSI  FUNCTION 
C  

DELT1= ( - BCB/2 . 0 ) * ( R+DD+LAM ( J , K ) *R ) 

DELT2=(-BCB/2.0)*(RR-EE-LAM(J,K)*RR) 

DELT3=( -P0/(2. 0*3. 1^1593 ))*(BETA/(LAM( J, K)**2-BETA**2)) 

DELTU= ( DELT1*PHI3+DELT2*PHI9+DELT3 ) * ( 1 . 0-FOU ) 

DELT5= ( DELT1+DELT2+DELT3 ) * ( 1 . 0+FOU ) 

DELT6=-PHIE(J,K)*CCE/CE 

IF(N.LE.U.O)  GO  TO  20U 

DELT7=CEXP(-PHI1) 

DELT8=CEXP(-PHI2) 

GO  TO  205 

204  DELT7=-CEXP(-PHI1) 
DELT8=-CEXP( -PHI2 ) 

205  DELT9=(PHlU*DELT7+PHI5*DELT8+PHl6)*(l .0-FOU) 
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DELTll=(PHll++PHI5+PHl6)*(l.O+FOU) 
DELT12= (DELT9+DELT11 ) * ( CCE/CE ) 

DELT(J,K)=(DELTU+DELT5)*((PEN-EIC(J,K))/CE)+DELT6+DELT12 
201  CONTINUE 

RETURN 
END 


$IBFTC  SUB3 

SUBROUTINE   STRESS(ll) 

C  THE  CALCULATION  OF  THE  STRESSES 

C  — - 

COMPLEX  DLJ,DLK,LAM(15,20),EI(15,20),EIC(15,20),P0,PHTE(15,20),DEL 
IT (15, '20) 
REAL  ZZ(U) 

COMMON/TRA/NN , MM , DLJ , DLK , N , C , BETA , CC , LAM , EI , EIC , PO , PR , DEJ , DEK 
C  QMMON/FUN/PHTE , DELT 

C0MM0N/STR/A(15,20).B(15,20),S(15,20),BMAX(15,20),AMAX(15,20),SMIN 
1(15,20),ALFA1(15,20),ALFA2(15,20) 
COMMON/POIS/ZZ 

F0U=ZZ(II)/(1.0-ZZ(II)) 

DO  301  J«=2,NN 

DO  301  K=1,MM 

IF(AIMAG(EI(J,K)).LE.0.0)  GO  TO  302 

A(J,K)=2.0*REAL(PHIE(J.K))+REAL(DELT(J,K)) 

B(J,K)=AIMAG(DELT( J,K) ) 

5(J,K)=2.0*REAL(PHIE(J,K))-REAL(DELT(J,K)) 

CAR=(A(J,K)+S(j,K))/2.0 

CHIP=(A(j,K)-S(j,K))/2.0 

BMAX(  J,  K)=3QRT  (CHIP**2-t-B(  J,  K)**2 ) 

AMAX ( J , K ) =CAR+BMAX ( J , K ) 
SMIN(J,K)=CAR-BMAX(J,K) 
GA=B(J,K)/CHIP 

ALFAl(j,K)  =  (ATAN(GA)/2.0)*(l80.0/3.1i+1593)+90.0 
ALFA2(j,K)=ALFAl(j,K)-90.0 
GO  TO  301 
302  A(J,K)=2.0*REAL(PHIE(J  K))+REAL(DELT(J,K) )+AIMAG(El(j,K)) 
B(J,K)=AIMAG(DELT(J,K)) 

S(J,K)=2.0*REAL(PHIE(J,K))-REAL(DELT(J,K))+AIMAG(EI(J,K))*F0U 
CAR=(A(j,K)+S(j,K))/2.0 
CHIP=(A(j,K)-S(j,K))/2.0 
BMAX(J,K)=SQRT(CHIP**2+B(J,K)**2) 
AMAX(J,K)=CAR+BMAX(J,K) 
SMIN(J,K)=CAR-BMAX(J,K) 
GA=B(J,K)/CHTP 
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ALFAl(j,K)=(ATAN(GA)/2.0)*(l8O.0/3-l!+1593)+90.0 
ALFA2( J, K) =ALFA1( J , K) -90 . 0 
301  CONTINUE 
RETURN 
END 

$IBFTC  SUBU 

SUBROUTINE  PHl(ll) 

C     THE  INTEGRAL  OF  PRTE  FUNCTION 
C     — 

COMPLEX  DLJ,DLK,LAM(15,20),CAM(15,20),EI(15,20),EIC(15,20),P0,PHT1 
1,PHT2,PHI3,PHI9,BCB,SPHI1, SPHI2, SPHI3,SPHIl+, SPHT5, SPHI6, SPHI7, SPHI 
28, SPHI9, SPHIIO, SPHI11 , SPHT12, SPHI13, SPHI1U , SPHI15 , SPHIA, SPHIB, SPHI 
321, SPHI22, SPHI23, SPHI2U , SPHT25, SPHI26, SPHT27, SPHI28, SPHI29, SPHI30, 
USPHI31, SPHI 32, SPHI33, SPHI3U, SPHI35, SPHIC , SPHE), SPHI(15, 20) , SHl(l5, 
520) 

REAL  N,ZZ(1+) 

COMM0N/TRA/NN,MM,DLJ,DLK,N,C,BETA,CC,LAM,EI,EIC,P0,PR,DEJ,DEK 

COMMON/PH/SPHI, SHI 

COMMON/POIS/ZZ 

F0U=ZZ(ll)/(l.0-ZZ(ll)) 

IF(N.LE.l+.0)  GO  TO  1+02 

PHIl=2.0*3.lUl593*(l-0/N)*PO 

PHI2=2.0*3.H+1593*(-1.0/N)*PO 

PKI3=CEXP(PHI1) 

PHI9=CEXP(PHI2) 

GO  TO  1+03 
1+02  PHI1=2.0*3.11+1593*(-1.0/N)*PO 

PHI2=2.0*3.LUl593*(l.O/N)*PO 

PHI3=-CEXP(PHI1) 

PHI9=-CEXP(PHI2) 
403   BC=1.0/((2.0*3.lUl593)*(l-0-BETA)) 

BCB=BC*PO 

DO  1+01  J=1,NN 

DO  1+01  K=1,MM 

IF(J.EQ.I.AND.K.NE.I)   GO  TO  1+01 

SPHI1=LAM ( J , K ) +BETA 

SPHI2=LAM(J,K)+1.0 

SPHI3=LAM(J,K)-1.0 

3PHlU=LAM( J , K) -BETA 

SPHI5=SPHI1*CL0G(SPHI1) 

SPHl6=SPHI2*CL0G ( SPHI2 ) 

SPHI7=3PHI3*CL0G ( SPHI3 ) 

SPHl8=SPHll+*CL0G(  SPHIU ) 

SPHI9=SPHI5-SPHl6+l .O-BETA 

SPHI10=SPHI7-SPHI8+1 . O-BETA 
.  SPHIU=(SPHIl/2.0)*SPHI5-(SPHIl**2)/l+.0-BETA*SPHI5+BETA*SPHIl-(SPH 
112/2 . 0 )  *SPHl6+ ( SPHI2**2 )/!+.  0+SPHI6-SPHI2 
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SPHI12=( SPHI3/2 .0 )*SPliI7- ( S?tiI3**2)/k .0+SPHI7-SPHI3-( SPHlU/2 . o)*SP 
lHl8+(SPHll+**2)/l+ .0-BETA*SPHl8-fBETA*SPHll+ 
SPHI13= ( -BCB/2 . 0 ) * ( SPHI9+SPHIll+( i . 0-BETA ) *LAM( J , K ) ) 
SPHI1U=(-BCB/2.0)*(SPHI10-3PHI12+(BETA-1.0)*LAM(J.K)) 
SPHI15=(-P0/(1+. 0*3. 1^1593) )*(SPHI8-SPHI5+2.0*BETA) 
SPHIA= ( SPHI13+SPHI1U+SPHI15 ) * ( 1 . 0+FOU) 
SPHIB=  ( SPKI13*PHI3+SPHIlit*PHI9+SPHI15  )*(  1 .0-FOU) 

cam(j,k)=real(lam(j,k))-aimag(iam(j,k))*po 
sphi21=cam( j , k) +beta 

SPRT22=CAM(J,K)+1.0 

SPHI23=CAM(J,K)-1.0 

SPHI2U=CAM( J , K ) -BETA 

SPHI25=SPHI21*CL0G ( SPHI21 ) 

SPHI26=SPHI22*CL0G( SPHI22 ) 

SPHI27=SPHI23*CL0G( SPHI23) 

SPHI28=SPHI2U*CL0G( SPHI2U ) 

SPHI29=SPHI25-SPHI26+1 .0-BETA 

SPHI30=SPHI27-SPHI28+1 . 0-BETA 

SPHI31=(SPHI2l/2.0)*SPHI25-(SPHI21**2)/l*.0-Bm*SPHI25+BETA*SPHI21 
l-(SPHI22/2.0)*SPHI26+(  SPHI22**2)/i+ .0+SPHI26-SPHI22 

SPHI32=( SPHI23/2 .0)*SPHI27- ( SPHI23**2 )/U .0+SPHI27-SPHI23- ( SPHI2U/2 
1 . 0 ) *SPHI28+( SPHI2U**2 ) /k . 0-BETA*SPHI28+BETA*SPHI2l+ 

SPHI33=(-BCB/2.0)*(3PHI29+SPHI31+(1.0-BETA)*CAM(J,K)) 

SPHI31+=(-BCB/2.0)*(SPHI30-3PHI32+(BETA-1.0)*CAM(j,K)) 

SPHI35=(-PO/(U.O*3 .1^1593) )*( 3PHI28-SPHI25+2 .0*BETA) 

SPHIC= ( SPHI33+SPHI31++SPHI35  )  *  ( 1 . 0+F0U ) 

SPHID=  ( SPHI33*PHI3+SP}n:  3i+*PHI9+SPHI  35  )  *  ( 1  •  0-FOU) 

SHI(J,K)=SPHIC+SPHID 

SPHI ( J , K)=3PHIA+SPHIB 

IF(J.EQ.l)   GO  TO  UOl 

SHI(J,K)=SHI(J,K)-SHI(1,1) 

SPHI(J,K)=SPHI(J,K)-3PHI(1,1) 
1+01  CONTINUE 

RETURN 

END 

$IBFTC  SUB5 

SUBROUTINE  PSl(ll) 

C     THE  MODIFICATION  OF  PHI  FUNCTION  FOR  THE  DISPLACEMENT 
C     

COMPLEX  DLJ,DLK,LAM(15,20),EI(15,20),EIC(15,20),PO,GI(15,20),ADCD, 
1CAM( 15 , 20 ) , ACD , ACE , PEN , BCB , ADD , AEE , PHI1 , PKT2 , PHI3 , PHI9 , APHlU , APHI5 
2 , APHI6 , APHI7 , APHI8 , APHIE 

REAL  N,ZZ(1+) 

COMM0N/TRA/NN,MM,DLJ,DLK,N,C,BETA,CC,LAM,EI,EIC,P0,PR,DEJ,DEK 

COMMON/GES/GI 

COMMON/POIS/ZZ 


FOU=ZZ(ll)/(l.O-ZZ(ll)) 
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IF(N.LE.U.O)   GO  TO  502 

PHT1=2 .0*3  .lifl593*(  1 .0/N)  *P0 

PHI2=2.O*3.lUl593*(-1.0/N)*P0 

PHI3=CEXP(PHI1^ 

PHI9=CEXP(PHI2) 

GO  TO  503 

502  PHIl=2.0*3.11tl593*(-1.0/N)*P0 
PKE2=2 .0*3 .1U1593*(1 .0/N) *P0 
PKE>-CEXP(PHIl) 
PHI9=-CEXP(PHI2) 

503  BC=1.0/((2.0*3.11+1593)*(l.O-BETA)) 
BCB=BC*PO 

DO  501  J=1,NN 
DO  501  K=1,MM 

if(j.eq.i.and.k.ne.i)  go  to  501 

cam(j,k)=real(lam(j,k))-aimag(lam(j,k))*po 

acd=(cam(j,k)**2-1.0)/(cam(j,k)**2-beta**2) 

ADCD=CL0G(ACD)/N 

ace=(1.0/(cc*c) )*cexp(adcd) 

fen=real(ei(j,k))-aimag(ei(j,k))*po 

add=clog((cam(j,k)+beta)/(cam(j,k)+i.o)) 

aee=clog((cam(j,k)-i.o)/(cam(j,k)-beta)) 

aphlu=(-bcb/2.0)*(add+cam(j,k)*add-beta+1.0) 

aphi5=(-bcb/2.0)*(aee-cam(j,k)*aee+beta-1.0) 

APHl6=(-PO/(i+.0*3.11+1593))*CLOG((CAM(j,K)-BETA)/(CAM(J,K)+BETA)) 
APHI7=(  APHlU*Pin;3+APHI5*PHI9+APHl6  )*(  1 .0-FOU) 
APHI8= ( APKCU+APHI5+APHI6 ) * ( 1 .0+FOU) 
APHIE= ( APHI7+APHI8 ) /ACE 

GI(J,K)=(REAL(APHIE)-AI>1AG(APHIE)*P0)*(EI(J,K)-PEN) 
IF(J.EQ.l)   GO  TO  501 
GI(J,K)=GI(J,K)-GI(1,1) 
501  CONTINUE 
RETURN 
END 

$I£FTC   3UB6 

SUBROUTINE  DISPLA(ll) 

C            THE  ELASTIC  DISPLACEMENT 
C  

COMPLEX  SPHI(15,20),SHI(15,20),GI(15,20),DLJ,DLK,LAM(15,20),EI(15, 
120),EIC(15,20),P0,DI 
REAL  N,ZZ(U) 

CCMMON/TRA/Nr^M-IjDLJjDLKjNjCjBETAjCCjLAMjEI^IC.PO^R.DEJjDEK 

COMMON/PH/SPHI ,  SHI 

COMMON/GES/GI 

CCMM0N/DIS/U(15,20),V(15,20) 

CCMMON/POIS/ZZ 

COMMON/CHIF/ADJ 
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F0U=ZZ(II)/(1.0-ZZ(II)) 

DO  601  J=2,NN 

DO  601  K=1,MM 

IF(AIMAG(EI(J,K)).LE.O.O)  GO  TO  602 

DI=((3.0-1+.0*ZZ(II))*SPHI(J,K)+SHI(J,K)-GI(J,K))*(1.0+ZZ(II)) 

GO  TO  603 

602  DI=((3.0-4.0*ZZ(II))*SPHI(J,K)+SHI(J,K)-GI(J,K)+PO*(1.0-FOU)*AIMAG 
1(EI(J.K))**2/2.0)*(1.0+ZZ(II)) 

603  U(J,K)=REAL(DI) 

V ( J , K ) =ADJ-AIMAG ( DI ) 
601  CONTINUE 
RETURN 
EI© 

$IBFTC  SUB7 

SUBROUTINE  VISCO(ll) 

C     THE  VISCOELASTIC  DISPLACEMENT 
C     

COMPLEX  DLJ.DLK,LAM(15,20),EI(15?20),EIC(15,20),P0 
REAL  N,ZZ(i+),AT(l8),VT(l8),UV(l8),W(l8) 

COMMON/TRA/M,MM,DI^DLK,N,C,BETA,CC,LAM,EI,EIC,PO,PR,DEJ,DEK 

C0MM0N/DIS/U(15,20),V(15,20) 

C  OMMON/VET/NL , DT , DET 

C        XXMUKXXKXXXXXMXMXHHKXXXKXMXXXMXKXXXXMKXXXKXKMXXHKXXIIKXXXXXMXX 

C  *  UV(LK)=  THE  VISCOELASTIC  DISPLACEMENT  in  horiz.  direction  * 

c  * * 

C  *  W(LK)=THE  VISCOELASTIC  DISPLACEMENT  IN  VERTICAL  DIRECTION* 

C  * * 

C  *  AT(LK)=  TIME  IN  HOURS  * 

C  * * 

C  *  VT(LK)=  TIME  IN  DAYS  * 

C  * * 

C  *  EleER/E  * 

C        XXXXttXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKKXXXXXXXXXXXKXXXX 

DO  701  J=2,NN 

DO  701  K=1,MM 

DO  705  LK=1,NL 

P=(LK+l)/2-l 

VeLK/2-1 

IF((LK-l)/2.EQ.LK/2)  GO  TO  702 

AT(LK)=DT*10.0**V 

VT(LK)=AT(LK)/2U.O 

GO  TO  703 

702  AT(LK)=DET*10.0**P 
VT(LK)=AT(LK)/2U.0 

703  El=8358.0/( (11112. 0)*(AT(LK)**0. 025953)) 
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uv(lk)=u(j,k)/e1 
w(lk)=v(j,k)/ei 
705  continue 

weite(6,27)  ei(j,k) 

write(6,28)  (at(lk)  ,vt(lk)  ,uv(lk)  ,w(lk)  .lk=1,nl) 

27  FORMAT(///lHO,UHX/H=,F9.5,7X,UHy/H=,F9.5) 

28  FORMAT (1H0,10H  TIME(HRS)  ?3X,10HTIME(DAYS)  ,6X,5HUV(X)  ,5X,5HW(Y)/// 
1(1X>E9.2,F12.5,3X,2F10.5)) 

701  CONTINUE 
RETURN 
END 


$DATA 
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APPEKDEC  D 

Influence  Diagrams 

This  Appendix  contains  influence  diagrams  for  vertical  normal  stress 
and  horizontal  normal  stress  along  vertical  sections  for  \i  =   0.3,  Ci  -   15°, 
30°,  k5° ,  60°  and  75*,  and  various  embankment  shapes.  Influence  diagrams 
for  horizontal  and  vertical  shear  stress  for  a  =  ^5°,  H  =  0.3  and  various 
embankment  shapes  are  also  given. 

The  diagrams  indicate  the  stress  due  to  the  embankment  weight  alone, 
stresses  due  to  the  weight  of  material  underlying  the  embankment  must  be 
superimposed  to  obtain  the  total  stress.  The  stresses  are  expressed  in 
dimensionless  form  as  a   /7H,  a   //"H,  or  *     /7H.  The  coordinates  are  also 
in  dimensionless  form.  For  convenience  in  the  semi-logarithmic  plot, 
the  depth  is  measured  from  the  top  of  the  embankment,  and  designated 
y/H.  This  is  in  contrast  to  the  discussion  in  the  body  of  the  paper 
where  the  vertical  distances  are  measured  from  the  base  of  the  embank- 
ment and  designated  y/H. 

Bach  of  the  four  diagrams  in  a  given  figure  refers  to  a  particular 
vertical  section,  shown  schematically  on  the  diagram.  The  upper  left 
diagram  indicates  stresses  along  the  center  line;  the  upper  right  dia- 
gram indicates  stresses  along  a  vertical  section  midway  between  the 
center  line  and  the  toe  of  the  slope;  the  lower  left  diagram  indicates 
stresses  along  a  vertical  section  through  the  toe  of  the  slope;  the 
lower  right  diagram  indicates  stresses  along  a  vertical  section  located 
a  distance  from  center  line  equal  to  1.5  times  the  distance  from  the 
center  line  to  the  toe  of  the  slope. 
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George  Y.  Baladi  was  born  in  Aleppo,  Syria  on  February  8,  1937  and 
is  presently  a  citizen  of  the  Syrian  Arab  Republic.  He  attended  Aleppo 
University  from  September,  1957  through  June,  19&1,  receiving  a  Bachelor 
degree  in  Civil  Engineering.  He  then  taught  soil  mechanics,  strength  cf 
materials,  and  mathematics  for  two  years  at  Aleppo  University. 

In  September,  I963,  he  started  his  graduate  work  at  Purdue  Univer- 
sity, and  graduated  in  June,  1965  with  a  Master  of  Science  in  Civil  En- 
gineering, majoring  in  soil  mechanics.  For  the  next  two  years  and  a 
half,  he  worked  on  his  Ph.D.  at  Purdue  University. 

Mr.  Baladi  is  an  Associate  Member  of  the  American  Society  of  Civil 
Engineers,  and  a  member  of  the  Syrian  Society  of  Civil  Engineers. 

In  1963,   Mr.  Baladi  was  married  to  Layla  Zeito  and  they  are  the 
parents  of  Zena  and  Joseph. 
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